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ABSTRACT 

We compute the absorption efficiency (Qabs) of forsterite using the discrete dipole 
approximation (DDA) in order to identify and describe what characteristics of crystal 
grain shape and size are important to the shape, peak location, and relative strength 
of spectral features in the 8-40 fim wavelength range. Using the DDSCAT code, we 
compute Qahs for non-spherical polyhedral grain shapes with a^g = 0.1 /nm. The shape 
characteristics identified are: 1) elongation/reduction along one of three crystallographic 
axes; 2) asymmetry, such that all three crystallographic axes are of different lengths; 
and 3) the presence of crystalline faces that are not parallel to a specific crystallographic 
axis, e.g., non-rectangular prisms and (di)pyramids. Elongation/reduction dominates 
the locations and shapes of spectral features near 10, 11, 16, 23.5, 27, and 33.5 /xm, 
while asymmetry and tips are secondary shape effects. Increasing grain sizes (0.1 - 
1.0 fim) shifts the 10, 11 /um features systematically towards longer wavelengths and 
relative to the 11 /xm feature increases the strengths and slightly broadens the longer 
wavelength features. Seven spectral shape classes are established for crystallographic a-, 
b-, and c-axes and include columnar and platelet shapes plus non-elongated or equant 
grain shapes. The spectral shape classes and the effects of grain size have practical 
application in identifying or excluding columnar, platelet or equant forsterite grain 
shapes in astrophysical environs. Identification of the shape characteristics of forsterite 
from 8-40 /Ltm spectra provides a potential means to probe the temperatures at which 
forsterite formed. 
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1. INTRODUCTION 



The silicate crystal mass fraction, /cryst, which is the mass ratio of silicate crystals to all 
silicates, is an analytical tool that can be used to study the taxonomy of comets and the evolution 
and radial stratification of proto-planetary disks. Crystalline silicates are rare in the ISM, < 2.2% 
along the line-of-sight to the Galactic Center and < 5% along other lines-of-sight (Li Sz Draine 



2001). In contrast, primitive Solar System materials reveal significantly higher /cryst- Thermal 
models of comet comae produce values of f cryst as high as ^ 0.8 (Harker et al. 2002, 2004) and 
generally comets have /cryst ^ 0.2 ( Kelley fc Wooden|2009 ). Stardust samples of comet 81P/Wild 2 



include large (~ 20 ^m) single mineral crystals (Brownlee et al. 2006 ) and X-ray microprobe analysis 



indicates a /cryst > 0.5 ( Ogliore et al.|200"9 ). These crystals are thought to be condensates ( |Brownlee| 
et al.|20()6 ; Nakamura et al.|2008t Tsuchiyama et al.|2009 ) . The discrepancy between ISM and Solar 
System abundances of crystalline silicates indicates the Solar System crystals must be a product 
of the formation and evolution of a dusty disk with a radial thermal gradient. The hot inner 
regions of proto-planetary disks have temperatures (T > 1000 — 1400 K) appropriate for rapid, 
direct gas phase condensation or annealing of amorphous ISM silicates into crystalline silicates 
( Bell et al.|[2000l |Gail||2002[ |2010t |Henning||2003t [Wooden et aL]|2007t |Henning||2010D . The presence 
of these hot refractory materials in cold, thermally unaltered comets indicates an efficient process 



of transporting materials from the hot inner regions of the disk to the cold outer regions ( Ciesla & 



Cuzzi 2007 Hughes Sz Armitage 2010). Further, the abundance of crystalline silicates at different 
heliocentric distances serves as a tracer of the efficiency of radial transport. Hence, /cryst provides 



constraints for radial transport in protoplanetary disk models ( Wooden et al.|[2007 Wooden]|2008 ) . 
In turn, radial transport in disks is used to interpret where, how, and at what time(s) and over 
what duration(s) the crystals formed and over what subsequent timescales they were transported 



outward into the comet forming zones of the Solar proto-planetary disk (Gail 2004 Ciesla 2010 



Hughes & Armitage 2010). 



The increasing sophistication of remote sensing observational techniques and advances in the 
analytic interpretation of infrared (IR) spectral energy distributions (SEDs) emitted by dust in 
comet comae, as well as constraints imposed by collected materials in comet sample-return missions, 
now necessitate modernization of astrophysical dust grain models to accurately derive /cryst- As the 
signal-to-noise quality of observations has improved over the years through space-based facilities 
such as the Spitzer Space Telescope, new subtleties in the IR silicate emission features evident 
in SEDs for cometary comae are discernible. Current modeling techniques are insufficient and 
cannot replicate many of the complex spectral signatures in the observed SEDs. Recent laboratory 
and modeling efforts suggest that to account for crystalline resonances, it does not suffice to only 
consider grain composition, but it is required to also consider grain shape and size. Grain shape 
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and size can change the wavelength location, shape, and relative strength of the crystalline silicate 



features present from 8-40 //m 


Fabian et al.|2001 


Min et al.|2005b 


Andersen et al.|2006 


Tamanai 


et al.||2006a|b 


Imai et al.||2009 


Mutschke et al.||2009 


Tamanai et al.||2009 


Koike et al.||2010 


). The 



magnitude of the changes in peak wavelength location due to the variation of grain shape in these 
laboratory experiments is on the order of the changes in peak location arising from the addition 



of up to 25% iron to an Mg-rich olivine crystal (Koike et al.||2003 Kelley &: Wooden 2009). Thus 



the grain shape of the crystals must be considered in parallel to external grain composition when 
modeling the dust in comet comae, especially when the derivation of /cryst is sought. The shape 
of the crystals in comet dust could be indicative of the process of crystal formation either by 



gas phase condensation (Tsuchiyama 1998 Bradley et al. 1983 Kobatake et al. 2008), selective 



evaporation (Tsuchiyama & Tachibana 1998 Takigawa et al. 2009), or annealing (Koike et al. 



2006 



2010). Furthermore, if variations in /cryst are used to interpret radial transport efficiencies 
or dust evolution processes in proto-planetary disks (PPDs) (e.g., Manoj et al. 2011 Hughes & 



Armitage 2010), accurate model descriptions of the dust, including the crystalline component, are 



required ( [Olofsson et"aL]|2009[ [Oliveira et al.||2011 ). 



In this paper, we discuss results of new computational modeling of polyhedral shaped forsterite 
(Mg2Si04, the magnesium-rich end-member of the olivine group) using the discrete dipole approx- 
imation (DDA). We investigate the characteristics of grain shape and grain size that produce 
subtleties in synthetic thermal IR SEDs. Such crystalline silicate features are now routinely ob- 
served by remote sensing techniques in comets and proto-planeteary disks conducted at moderate 
spectral resolutions {R = AA/A > 1000) and high sensitivity (continuum flux densities of mJy to 
fJ-Jy). We present the initial stages of building a bottom- up model that will compute absorption 
efficiencies for large, porous, multimineralic aggregates with crystalline inclusions for the 8-40 /im 
IR spectral range. We systematically vary the axial ratios and external shapes to enable an in depth 
study of what characteristics of grain shape are important to the crystalline resonant features in 



the absorption efficiency spectra for forsterite, which is a dominant mineral in comet dust (Harker 



et al. 2007 


Lisse et al. 


2007 


Reach et al. 


2010 


Harker et al. 


2011 



crystalline resonance spectral features in the 8-40 region. The crystal shapes selected from this 
study provide realistic model crystals that can be incorporated in the future into a large, porous, 
multimineralic aggregates used in thermal dust models. 

A major challenge in deriving the silicate mineralogy of comets is understanding how the 
anisotropic nature of forsterite crystals affects the 8-40 ^m spectral features' shape, strength, and 



peak position. As forsterite is the main crystalline silicate identified in cometary comae (Campins 



Ryan|[T989l [Crovisier et"aLl[T997l [Wooden et ar][T999l [Harker et alT][2002l [20041 [Lisse et al. [[20061 



2007), accurate models for forsterite's emission/absorption features are the primary requirement 
to determine the silicate mineralogy of comets. However, there exists a dearth of models that fully 
account for the anisotropic nature of forsterite. Here we present a new modeling investigation into 
the 8-40 /im absorption resonance features of forsterite using the discrete dipole approximation 



(DDA) code DDSCAT (Draine fc Flatau 1994 Draine k Flatau 2008a), which fully accounts for 
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the anisotropic nature of forsterite. 

The paper proceeds in the following manner. In § [2| we briefly summarize three commonly 
used computational approaches to modeling grains and explain why DDA is the method chosen 
to compute polyhedral shaped forsterite crystals. In § [Sj we describe the details of how DDA is 
applied to our study of grain shape characteristics, and what constraints apply to using the publicly 
available DDA code DDSCAT ( |Draine Flatau||1994[ [Draine k Flatau||2008a| ). Section |4] presents 
the DDSCAT-calculated absorption efficiencies (Qabs)- We compute Qabs for three grains shape 
characteristic exercises: 1) the elongation/reduction of a single crystallographic axis; 2) rectangular 
prism shapes that have three different side lengths, and are therefore asymmetric; and 3) shapes that 
have different external geometries but have the same crystallographic axis ratio with an emphasis 
on exploring the effects of (di) pyramidal structures through a comparison of Qabs between prism - 
dipyramidal pairs. Also in § |4] we compute Qabs for variations in grain size from 0.1 - 3.0 //m in 
effective radius. In § [Sj we expand upon the results of § |4] and discuss the implications of crystal 
shape characteristics and their potential to act as a probe of formation enivrons, and we provide 
a practical application to the results to the case of cometary comae. In § |6] we summarize the 
major conclusions. The major forsterite bands we investigate are those located near 10, 10.5, 11, 
16, 19, 23, 27, and 33 /xm. In addition, we examine more subtle emission features appearing at 
approximately 11.9 and 25.2 /xm. 



COMPUTATIONAL APPROACHES TO OPTICAL CONSTANTS 



There are three broad approaches in modeling methods to compute the optical properties 
arising from an ensemble of small (less than a few microns in size) discrete grains that are sub- 
sequently incorporated into thermal grain models to replicate observed SEDs (for a review see 
Min 2009). The first approach is to assume that the grains are homogeneous in composition and 



spherical in shape. In this case, Mie theory (Mie 1908) can be used, and computations are simple 



and fast. However, in cometary samples from interplanetary dust particles (Molster et al. 2002 



Bradley||2003 e.g., ) crystalline grains are non-spherical. Forsterite grains in Stardust samples can 



have euhedral, polyhedral with well-defined sharp faces and edges, shapes ( Nakamura-Messenger 



eraI][20T2| ) or irregular shapes (i.e., not even or balanced in shape or arrangement) where irregular 



shapes may arise due to capture in aerogel (Zolensky et al. 2008). Yanamandra-Fisher & Hanner 



(1999) and Fabian et al. (2001) demonstrate that Mie computations using spherical particles in- 



adequately replicate spectrophotometric observations if crystalline silicates are present within the 
emitting dust ensemble. The symmetry of a sphere creates resonance effects at particular values of 



the index of refraction that are not seen in natural grains with irregular shape ( Yanamandra-Fisher 



k Hanner 1999 Fabian et al.||2001 Min 2009 Henning k Mutschke 2010). 



The second modeling methodology is a statistical approach to account for irregular, non- 
polyhedral, grain shapes. In this approach, an arbitrarily shaped particle is represented by a 
distribution of simple particles that have exact solutions. For example, the continuous distribution 
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of ellipsoids (CDE) method ( [Fabian et al]|2001[ |Min et al.||2003t [Bohren fc Huffman| [19831 ) is 
frequently used. This adopts a statistical distribution of ellipsoids, which may include spheres, that 
can vary in size (within the Rayleigh limit), principle axes ratio (or for crystals, crystallographic 
or optical axes), and orientation. The absorption characteristics for each ellipsoid can be exactly 
computed, and the grain shape irregularity is defined by the choice of ellipsoids. Within the 
statistical approach, ellipsoids are but one of several simple shapes with an exact solution that can 
be chosen. Other notable simple shapes that can be used to acquire the absorption characteristics 



of an irregularly shaped grain are hollow spheres (Min et al. 2003, 2005a) and gaussian random 



spheres (Muinonen et al. 1996: Volten et al. 2001 Mutschke et al. 2009). 



The third modeling methodology seeks to directly calculate the optical properties for arbitrarily 



shaped grains. These computationally intensive approaches are the DDA (Purcell &i Pennypacker 



1973) method, the coupled dipole approximation (CDA) ( |Min et al.|2008 ), or the T-Matrix method 



(TMM) ( [Mishchenko et al.[[l996[ ). The TMM ( [Mishchenko et al.[[2007| ) expands the incident and 
scattered electric fields in series of suitable vector spherical wave functions, wherein the relation 
between the columns of the respective expansion coefficients is established by means of a transition 
matrix (or T- matrix). This approach can be applied to the entire scatterer as well as to separate 
parts of a composite scatterer ( Mishchenko et al.|2007 ). The TMM has been successful in modeling 
the polarization of light scattered by aggregate grains ( Kolokolova et al.|2007 Kolokolova &; Kimura 



2010). In theory the TMM can account for anisotropic materials (Schmidt Sz Wriedt 2009), however 



to date the TMM has not been extended to model a discrete tri-axial forsterite crystal due to a 
reliance upon a reformulation of the TMM equations that are specific to a particular grain shape. 



The DDA (Purcell & Pennypacker 1973 Draine & Flatau 1994 Min et al. 2006) is an ab- 



sorption and scattering modeling method that allows for grains of arbitrary shape and crystalline 
composition such that all three crystallographic axes are accounted for simultaneously. With DDA 
the user can specify the geometry of the target (model grain), the size of the target, the alignment 
of multiple indices of refraction with respect to the target geometry, the wavelength range over 
which the model is to run, and the composition of individual dipoles. With this high level of con- 
trol, realistic targets of discrete crystalline grains can be constructed, permitting calculation of the 
target's absorption and scattering efficiencies. 

Our goal in this paper is to identify and describe characteristics of grain shape that affect the 
absorption of forsterite crystals in the context of forsterite's anisotropic shape (i.e., three distinct 
indices of refraction associated with three crystallographic axes). This necessitates a careful control 
on the exact shape of the target and a method that allows for computations to be done incorporating 
all three crystallographic axes simultaneously. The DDA approach satisfies both of these needs and 
is therefore selected for our study. 
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METHODOLOGIES 



The DDA technique is used to model crystalhne forsterite (Mg2Si04) between 8-40 iim with 
a wavelength spacing of AA = 0.05 /im. The DDA method represents an arbitrarily shaped particle 
by a collection of dipoles arranged on a three-dimensional lattice. In essence, each dipole is a polar- 
izable element that is a coarse representation of a large number of actual atoms within the target, 
such that with the appropriate indices of refraction the bulk dielectric constant of the material is 



reproduced by the dipoles (Purcell & Pennypacker 1973). Each dipole reacts to incident radiation 



producing a dipole field that in turn then interacts with all other dipoles in the DDA lattice or 
grid of dipoles that represents of the grain. This method is computationally intensive, yet also the 
most accurate for calculating absorption and scattering properties of grains. We use the publicly 
available DDSCAT code (Draine Sz Flatau 2008b| Version 7.0), which has been parallelized to run 
on the NASA Advanced Supercomputing (NAS) facility, Pleiades. 



3.1. DDA Targets 

Since the user specifies the location of each individual dipole in DDSCAT, any target shape 
can be created. This allows us to create a variety of target shapes that can be used to explore three 
main shape characteristics: 1) the elongation or reduction of a crystal along a single crystallographic 
axis; 2) shapes that differ in length for all three crystallographic axes, (i.e., crystallographic axial 
ratio), or rather are tri-axially asymmetric; and 3) the presence of faces and edges that are not 
aligned with a single crystallographic axis. The targets used to explore these shape characteristics 
are polyhedral in shape and include simple prisms and more complex dipyramidal structures that 
can be combinations of prismatic shapes with pyramidal tips. 

In addition to user-specified dipole location, any index of refraction with up to three differing 
crystallographic axes can be applied on a dipole-by-dipole basis, allowing for a unique composition 
to be applied to each dipole. Forsterite is anisotropic such that it has three distinct crystallographic 
axes, a, /3, 7, commonly referred to as the a-, b-, and c-axes. When creating our DDA targets, 
we are able to assign each dipole three indices of refraction that are each associated with one 
of forsterite's crystallographic axes. The crystallographic axes in the targets are assigned such 
that they are parallel to the x-, y-, and z- geometric axes. Using a variety of assignments, we 
are able to create targets that are elongated or reduced along our choice of crystallographic axis. 
These assignments allow for the investigations into the spectral effects associated with precise 
crystallographic axis ratios. 

A selection of the types of DDA targets used in this study are provided in Fig. [TJ The 
entire sample of polyhedral shapes includes shapes elongated along a single crystallographic axis 
(columnar), reduced along a single crystallographic axis (platelet), with all three axes unequal in 
length (asymmetric), prismatic shapes, (di)pryamidal shapes, and elongated dipyramidal shapes. 
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Elong ated Triangular 
Dipyramid (ETDP) 




a : b : c = 
0.47 : 1 . 0.59 
Brick 




Triangular 
I3ipyramid (TDP) 




TetrahedrDn 



a : ts : c = 
0.47 : 1 : 0.59 
Dflsyramidal Brick 




Elongated 
Hexagonal 
D^ryamid (EH DP) 




Trfangular 
Prism (Tri. Prism) 




a : b : c = 
0.47 : 1 :0.59 
Ellipsoid 



OcEalnodron 
(Octa) 




c-Shortened 
Fors'efile Dimension 
Brick Platelet 



Fig. 1. — A selection of some of the target shapes used in this study. Each sphere indicates the 
placement of a dipole in DDSCAT. The listed axial ratios for the ellipsoid, brick, and dipyramidal 
brick are a : b : c = 0.47 : 1 : 0.59, which is the axial ratio for the unit cell of forsterite. 
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3.1.1. Numerical Constraints 

There are two main numerical constraints that apply to the DDSCAT code. The first is that 
the lattice spacing d between dipoles must be small enough, and the number of dipoles used large 
enough, to accurately describe the shape of the target. This constraint is key for this work because 
we use a variety of polyhedral grain shapes. As such, we use a sufficiently large number of dipoles 
so that the target has relatively smooth faces and well-defined edges. Each dipole in the lattice is 
represented by a cubical volume element with side length, d. Thus, a sufficient number of these 
cubical volume elements are used to ensure that each cubical volume element is small compared 
to the overall size and shape of the target. The second constraint is stringent, and dictates that 
in order for the DDA to apply, d must be small compared to the wavelength of the radiant energy 
passing through the specified medium characterized by a complex index of refraction, m. This is 
satisfied by the condition \m\kd < 1.0 for Qabs calculations, where k = 27r/A is the wavenumber 



(Draine Sz Flatau 1994). If scattering calculations are needed, the condition is tighter, such that 



\m\kd < 0.5. 

A more useful form of the latter constraint relates |m|/c(i < 1.0 to the number of dipoles that 
are required to sufficiently model a target grain of a certain grain size. DDSCAT models a target 
on a three-dimensional lattice, where each cubic volume element with side length d is assigned to 
be a dipole with a particular index of refraction. Thus, each dipole element has volume and the 
total target volume is y = A'dipfi^, where N^ip is the total number of dipoles in the target. This 
volume must be equivalent to the actual physical volume of the polyhedral target. If the volume 
of the polyhedral shaped target were collapsed into a sphere of equivalent volume, then the radius 
of that sphere defines the effective radius of the target, agg, so that Oeg = d [(3A''dip)/(47r)]-^/^. 
Expressing Oeg in terms of A'^jip and \m\kd, yields 



/3iVdip\^ A 

«efr = ^.TT^ i—r m A;(i 1 
\ 327r^ / \m\ 

where A is the wavelength of radiant energy being considered. Applying the requirement for accurate 
absorption efficiency calculations, \'m\kd < 1.0, gives us a relation for the number of dipoles required 
for a target as a function of the incident radiation's wavelength, the complex index of refraction 
for the target material, and for the grain effective radius: 



3 



(2) 



Computing Qabs for a minimum wavelength of 3 /xm and assuming an index of refraction for 
forsterite at 3 fim of 1.8851 for the real part of the index of refraction ( Steyer|["l974 ) implies a target 
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with an effective radius of 1 fj,m requires 258 dipoles; 2 /im requires 2,063 dipoles; 3 requires 
6,961 dipoles; and 5 /im requires 32,224 dipoles. For grain radii QcS ^ 3.0 fim, we utilize roughly 
2,000 - 10,000 dipoles, depending on the target, for the full 3-40 fim wavelength range. 



3.1.2. Forsterite Spectral Resonances 



We seek to delineate how varying the shape and size characteristics of forsterite crystals affects 
the spectral crystalline resonant features' shape, peak location, and relative strength. By Kirch- 
hoff 's law of thermal radiation, the spectral features present in the absorption efficiencies (Qabs) will 
be present as emission features when these same grains are in radiative equilibrium with sunlight 



Wooden et al. 


1999 


Harker et al. 


2002 


2004 



As discussed in § [T| Mg-rich crystalline silicates are the dominant crystalline silicate component 



seen in comet C/1995 01 (Hale-Bopp) ( Crovisier et al. 1997 Harker et al. 2002 Lisse et al. 2006) 



and in comets in general ( Hanner &: Zolensky||2010 ). Resonant peak locations have been shown by 



observations (Molster et al. 2002) and by laboratory experiments to vary depending on the shape 



and size of the sampled grains ( |Min et al.|[2005bt [Mutschke et al.||2009[ [Tamanai et"ar]|2006a[ |2009 



Koike et al.|2010[ ). As such, we focus our study on crystalline forsterite with polyhedral shapes, and 



choose optical constants from Steyer (1974) to model our grains. The optical constants of Steyer 
(1974) were derived from specular reflectance measurements of a polished crystal surface and so 
inherently do not contain grain shape-dependent effects. Over the 8 -40 /xm wavelength range. 



the optical constants of Steyer ( 1974 ) have been proven to be equivalent to the optical indices of 



refraction derived by Fabian et al. (2001) from IR transmission measurements of polycrystalline 
forsterite powders (smaller than 1 fim) prepared by ball- grinding and dispersed in KBr. Note that 
the 49 jUm and 69 /im features were discerned from the measurements of polycrystalline powders 
(Fabian et al.||200l), but these features are beyond the wavelength range of this study. 



4. RESULTS 
4.1. Crystal Shape 

We investigate the effects of forsterite crystal shape on both the resonant absorption features' 
shape and the resonant features' peak location. The primary features of interest are the strong 
features located near 10, 10.5, 11, 16, 19, 23, 27, and 33.5 /im and the weaker, yet significant 
features near 10.5, 11.9, and 25.2 /im. The investigation focuses upon understanding how the 
crystallographic anisotropic nature of forsterite affects the 8-40 //m crystalline resonance features 
with respect to crystallographic axis length and specific grain shape characteristics, such as the 
importance of the presence of and angles between the vertices, edges, and faces of polyhedral crystal 
shapes. We begin by characterizing the spectral trends associated with elongating or reducing one 
of the crystallographic a-, b-, or c-axes for rectangular prisms. The rectangular prisms (henceforth 
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referred to as bricks) are the simplest polyhedral shape to examine the effects of elongating or 
reducing a crystallographic axis for the DDA since they provide a uniform number of dipoles along 
all the crystallographic dimensions. The spectral effects of elongation and reduction of bricks thus 
provide a fiducial for how the balance between crystal shape and length of the crystallographic axes 
affect the forsterite absorption resonance features. 



In § 4.1.1 , the absorption feature effects of elongations or reductions of a single crystallographic 
axis are characterized using the base shape of a cube. The equal lengths, or symmetry, of all the 
dimensions of a cube limit the effects on the absorption features due to the other two symmetric 
axes. This allows for identification of spectral trends associated with elongating and reducing the 
a-, b-, and c- crystallographic axes. However, the equal length symmetry of the two unaltered axes 
may be its own source of spectral effects. To test whether this is the case and confirm the spectral 



trends, we repeat the elongation and reduction experiment of bricks in § 4.1.2 using a base shape 
of a triaxial brick, where the two unaltered axes are set to the ratio of 0.8 and 1.2 times the length 
of the axis to be altered. We find that asymmetry does induce significant spectral effects, such that 
the two asymmetries (0.8,1.2 and 1.2,0.8) have divergent spectral features. Regardless, the spectral 
trends demonstrated by the symmetric case still hold for each asymmetric case. 

The divergent spectral behavior demonstrated by the asymmetric bricks indicates that the spe- 
cific crystallographic axial ratio is important to the overall absorption characteristics of a forsterite 



grain. In § 4.1.3, non-brick shapes with similar axial ratios are compared to demonstrate that the 
effects of elongations and reductions hold for more complex grain shapes and to provide a method 
to test the effects of differences in grain shape. Specifically, the comparison allows for a direct way 
to investigate the effects of specific grain shape with respect to the presence and angles between 
faces, edges, and vertices through examination of the absorption feature differences between shapes 
with exact or very similar axial ratios. 



4.1.1. Variation of Crystallographic Axes Lengths: The Symmetric Case 

Through understanding the effect of the simplest shape variations, it is possible to understand 
more complex shape effects. Hence, we start with elongations and reductions of the symmetric 
case of the rectangular prism: the cube. Shapes delineated by faces parallel to crystallographic 
axes are the simplest to relate to elongations or reductions of the axes lengths of a cube. Shapes 
delineated by two (triangular prism) or more faces being non-parallel to the crystallographic axes 
show similar trends as the elongated cube in the wavelengths and strengths of the features but the 
feature shapes are not predictable as a linear combination of elongations of a cube. 

The symmetry of a cube provides for the most controlled examination of the spectral effects 
related to an elongated or reduced single crystallographic axis. The elongations and reductions 
are defined by systematically varying the length of one of the crystallographic axes (the a-, b-, or 
c-axis) with respect to the other two while keeping the size of the grain constant at Oefr = 0.1 /im. 



- 11 - 



The elongations used are 1.2, 1.4, 1.6, 1.8, 2.0, 2.6, 3.0, and 4.0 times the original length of the 
cube-shaped crystal. For example, a 4.0 times crystallographic a-axis grain is essentially four cubes 
stacked upon each other in the a-axis direction giving ana:b:c = 4:l:l axial ratio. The 
reductions used are 0.8, 0.6, 0.4, 0.2, and 0.1 times the original length of a the cube-shaped crystal, 
such that at large reductions - 0.2 and 0.1 times - the grain shape is best described as a platelet. 

The absorption efficiencies of the elongations and reductions are displayed in Figs. [2] and |3] for 
the 8-13 fim and 13 - 40 //m range, respectively. The fiducial spectrum of the cube is described 
as follows: the '10 ^m feature' and '11 ^m feature' appear at 9.8 and 11.0 /xm, respectively, with 
minor peaks at 10.2 and 10.6 ^m. The stronger 18 and 23 /im features are asymmetric, such that 
the 18 /im feature has two shoulders on the long wavelength side of the peak, and the 23 /im is 
a triple-peaked, 'trident '-shaped feature. The weaker 16, 27, and 33.5 /xm features are singular 
symmetric peaks. 

The elongations/reductions of any one axis exhibit consistent spectral trends, i.e., the 2.0 times 
elongations are intermediary between the 1.0 and 4.0 times elongations. The effects of elongations 
and reductions approach having distinct spectral features as the elongations and reductions become 
'extreme', i.e., > 2.6 times elongated or < 0.2 times reduced. In these 'extreme' cases, we refer to 
elongated shapes as columns (occasionally referred to as 'whiskers' in the literature) and reduced 
shapes as platelets. 

Below, we discuss the specific spectral effects with respect to feature shape, strength, and 
location for elongations/reductions for each of the a-, b-, and c-axes. 

Elongating the cube along the a-axis shifts the 10 /im peak to shorter wavelengths and the 11 nm 
peak to longer wavelengths (Fig. [2| left panel, top). The 10 fim. feature also slightly decreases in 
relative strength with respect to the 11 /Lim feature (hereafter, all references to changes in relative 
strength are relative to the 11 //m feature). For the weaker features in the 8-13 //m range, 
elongating the a-axis causes the 10.6 //m feature to become a distinct separate peak near 10.5 
that shifts to shorter wavelengths and decreases in feature strength; the 10.2 and 11.9 /im peaks 
become more distinct with increased relative strengths. The cube has a 23 /im 'trident'-shaped 
feature that with a-axis elongation morphs into a bifurcated pair of peaks with a weak third peak 
on the long wavelength end of the feature (Fig. [3| left panel, top). This bifurcated pair of peaks is 
unique to a-axis elongations, and should be discernible in spectroscopic observations. 

Elongating the cube along the b-axis shifts the 10 /xm peak to shorter wavelengths and the 
11 /im peak to longer wavelengths (Fig. [2| middle panel, top). For the weaker features in the 
8 -13 /im range, elongating the b-axis causes the 10.6 /xm feature to split into two narrow features 
at 10.8 /xm and at 11.0 /xm, with the 10.8 /xm increasing and the 11.0 /xm decreasing in relative 
strength. With increasing b-axis elongation, the 11.9 /xm feature becomes an increased absorptivity 
tail for the red side of the 11 /xm feature. The cube has an 18 /xm feature with two long wavelength 
shoulders that with b-axis elongation fade and a short wavelength-side shoulder develops, becoming 
a minor companion peak for 'extreme' elongations (Fig.lsl middle panel, top). The cube has a 23 /xm 
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Fig. 2. — The 8-13 /im absorption efficiencies normalized to the 11 fim peak for a- {left panel), 
b- {middle panel), and c-axis {right panel) elongations and reduction of the crystallographic axes 
with respect to a fiducial cube shape {black, middle). The top curves in warm colors are 1.2, 1.4, 
1.6, 1.8, 2.0, 2.6, 3.0, and 4.0 times the original length of the cube elongations. The bottom curves 
in cool colors are 0.8, 0.6, 0.4, 0.2, and 0.1 times the original length of the cube reductions. 
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Fig. 3. — The 13 - 40 /iin absorption efficiencies normalized to the 11 nm peak for a- {left panel), 
b- {middle panel), and c-axis {right panel) elongations and reduction of the crystallographic axes 
with respect to a fiducial cube shape {black, middle). The top curves in warm colors are 1.2, 1.4, 
1.6, 1.8, 2.0, 2.6, 3.0, and 4.0 times the original length of the cube elongations. The bottom curves 
in cool colors are 0.8, 0.6, 0.4, 0.2, and 0.1 times the original length of the cube reductions. The 
Qabs values are multiplied by a factor of 1.5 for visual clarity. 
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'trident '-shaped feature that with b-axis elongation becomes a well-defined single peak plus a minor 
peak between the 18 and 23 features. The ability to discern b-axis elongation from a- or c-axis 
elongations is possible through the absence of the distinct spectral features associated with the a- 
or c-axis elongations. 

Elongating the cube along the c-axis shifts the 10 iim feature to longer wavelengths for non- 
'extreme' elongations. For 'extreme' elongations, the 10 /im feature becomes a shelf feature that 
intersects with the 10.6 //m feature. As elongation increases, the 10.6 //mfeature shifts to 10.5 fim 
and becomes a distinct peak that is nearly equal in strength to the 11 fim feature (Fig. [2| right 
panel, top). The 11 ^m feature shifts to slightly shorter wavelengths, which is the opposite behavior 
of the a- and b-axis elongations. With c-axis elongation, the 16 fim feature significantly increases 
in relative strength and the 18 nm feature morphs from a single peak into a bifurcated, double 
peak. The cube has a 23 iim 'trident '-shaped feature that with c-axis elongation becomes a single 
prominant peak with increasing strength and shifts slightly toward longer wavelengths; by 4.0 times 
elongation the 23 fim feature is nearly as strong (approximately two-thirds) as the 11 /imfeature 
(Fig. |3| right panel, top). With c-axis elongation, the 33.5 fim. feature does not change in feature 
shape, location, or strength. The c-axis elongation behavior for the 16 and 33.5 /im features is 
the opposite for what is observed for a- and b-axis elongations. The 27 fim feature also increases 
in relative strength and develops a distinct 'shark-tooth' shape, rising quickly and falling slowly. 
The c-axis 'extreme' elongations produce the combination 10 ^m shelf - 10.6 /im - 11 /xm feature, 
significantly increased 23 micron feature strength, and a 'shark-tooth' 27 /xm feature, which are 
distinct spectral patterns for c-axis elongations compared to a- and b-axis elongations that should 
be discernible in observations. 

In contrast to elongations, where one axis becomes increasingly longer over the other two, 
the case of reductions is where one axis becomes shorter than the other two. As exhibited by the 
spectra at the bottom of Figs. [2] and [3j axial reductions with one axis significantly shorter than 
the other two do not follow the same trends established with elongations. In general, reductions 
create more complex spectral feature trends with increasing reduction, which can be considered a 
combination of elongations of the two other (unaltered) axes. 

Reducing the cube along the a-axis from 0.8 times to 0.1 times the original length of the 
cube shifts the 10 and 11 ^m features to longer wavelengths (Figs. [2] left panel, bottom). The 
10.2 /im feature increases in relative strength and merges with the 10 //m feature for 'extreme' 
reductions (< 0.2 times) creating a rounded peak. Also at 'extreme' reductions, the cube's 10.6 fim 
feature moves to shorter wavelengths and becomes a narrow peak at 10.5 /xm with increasing feature 
strength. The combination of the 10.2 and 10.6 //m features' shape and location is similar to that of 
c-axis elongation. However, at 'extreme' reductions, the 11 //m feature develops an absorption tail 
that is similar to b-axis elongations. The 16, 18, 23, and 33.5 fim features all shift to significantly 
longer wavelengths (AA ^1.0 /im) and become relatively stronger with increasing a-axis reduction; 
these feature shifts are characteristic of b- and c-axis elongations. A-axis reductions also create a 
significant increase in the relative strength of the 23 /xm feature, and a 'shark-tooth'-shaped 27 /im 
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feature, which are characteristic of c-axis elongations. The overall feature structure suggests that 
the spectral features for a-axis reductions are in reality a combination of b- and c-axis elongation 
spectral features. The pattern of axial reductions of one crystallographic axis having spectral 
characteristics similar to a combination of elongations of the two unaltered axes holds for a-, b- 
and c-axis reductions. 

Reducing the cube along the b-axis shifts the 10 and 11 jim features to longer wavelengths. Sim- 
ilar to a-axis elongations, b-axis reductions have a distinct 10.6 /xm feature for all non-'extreme' 
reductions (> 0.4 times). Towards 'extreme' b-axis reductions, the 10 and 10.6 /xm features be- 
come a single, rounded peak centered near 10.2 /xm. With b-axis reductions, the minor 11.9 /um 
feature becomes increasingly strong, which is a characteristic of a-axis elongations. With b-axis 
reductions, the 16, 18, 23, 27, and 33.5 fim features all shift to longer wavelengths. Also, the 18 
and 23 jim features significantly change in feature shape. The cube has a 18 //m single peak with 
two long wavelength shoulders that with b-axis reductions morphs into a bifurcated double-peak 
that decreases in relative strength. The cube has a 23 jum 'trident'-shaped feature that morphs to 
a singular peak with significantly increasing relative strength. The 18 /xm bifurcated double-peak 
and the strong 23 /xm singular peak are both characteristics of c-axis elongations. 

Reducing the cube along the c-axis shifts the 10 fim feature to shorter wavelengths and shifts 
the 11 iim. feature to longer wavelengths. C-axis reductions are the only axial reductions that shift 
the 10 fim feature to shorter wavelengths, which is opposite the characteristic of both a- and b-axis 
elongations where this feature shifts to shorter wavelengths. As the c-axis is reduced, the 10 /xm 
feature also significantly decreases in relative strength, and the 10.6 fim feature diminishes to a 
small peak at 11.0 fim. Neither of these spectral changes are apparent in a- or b-axis elongations 
suggesting that not all spectral characteristics associated with a single axis reduction can be un- 
derstood as a combination of elongations of the two unaltered axes. With c-axis reductions, the 
18, 23, 27, and 33.5 ^m features shift to longer wavelengths; the 16 fim feature shifts to shorter 
wavelengths and significantly diminishes in relative strength. The cube has an 18 /xm feature with 
two long wavelength shoulders that with c-axis reductions morphs into a singular peak with a minor 
sub-peak on the short wavelength-side. The cube has a 23 /xm 'trident'-shaped feature that morphs 
into to a stair-stepped-shaped triple peak of increasing strength at longer wavelengths. The change 
in the 18 una feature's shape is similar to b-axis elongations, but the change in the 23 /xm feature's 
shape is not apparent in either a- or b-axis elongations. 

The more complex behavior for crystallographic axis reductions largely can be understood as 
a combination of elongations of the two unaltered axes. As the target axis is reduced, the other two 
axes becoming significantly longer, or elongated, in comparison. These combinations create distinct 
spectral features for a-, b-, and c-axis elongations that should be observationally identifiable in the 
spectra of astrophysical objects. 

In this section, the spectral changes presented for the symmetric case of the cube demonstrate 
that elongations and reductions of each crystallographic axis have distinct spectral characteristics 
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over the 8-40 fim spectral range. The trends with increasing elongation and reduction are 
predictable and different for each crystallographic axis. To summarize: The a-axis elongations 
have the distinct 10.5 /im feature and a bifurcated 23 fim double-peak. The a-axis reductions 
trend toward a 10.2 /xm rounded peak paired with a narrow 10.5 /im peak, and with increasing 
a-axis reduction the 15 - 40 //m features all move to significantly longer wavelengths and the 23 fim 
feature strength becomes significantly enhanced. The b-axis elongations have spectral features that 
are singular peaks, exhibit a 10.8 /im shelf, and an absorptive tail on the long wavelength side of 
the 11 /im feature. The b-axis reductions trend toward a singular, broad 10.2 /xm feature, and an 
enhancement in the 23 /xm feature strength that is slightly weaker than that of a-axis reductions. 
The c-axis elongations trend toward a 10 /im shelf, a bifurcated 18 /xm feature, and a 23 /im feature 
that is significantly enhanced in relative strength. The c-axis reductions exhibit a 10 /im feature 
that diminishes in relative strength and a 'stair-stepped' 23 /im feature. In the case of 'extreme' 
elongations (> 2.6 times the original length of the cube, also called 'columns'), and reductions (< 
0.2 times the original length of the cube, also called 'platelets'), the spectral features are distinct 
enough such that they should be discernible in observations. 



4-1-2. Variation of Crystallographic Axes Lengths: The Asymmetric Case 

The previous section detailed the spectral effects due to the elongation/reduction of a sin- 
gle crystallographic axis with respect to the base shape of a cube. Here, we repeat the elonga- 
tion/reduction exercise with an asymmetric rectangular prism (denoted as an asymmetric brick) 
with each of the three edges having differing lengths. 

The shapes for this exercise of the elongations/reductions of the asymmetric case have edge 
lengths of 0.8/, I.OIX, and 1.2/, where / is the edge length of the cube (in dipoles) in the previous 
section, and X is the multiplicative factor for the elongations and reductions. Since the edge length 
of the cube, c, is in every term, we drop it from the notation. For the 'fiducial bricks' {X = 1), there 
are six asymmetric realizations that include: two a-axis bricks, one with a longer b-axis (denoted 
a : b : c = 1.0 : 1.2 : 0.8) and one with a longer c-axis (a : b : c = 1.0 : 0.8 : 1.2); two b-axis bricks 
(0.8 : 1.0 : 1.2 and 1.2 : 1.0 : 0.8); and two c-axis bricks (0.8 : 1.2 : 1.0 and 1.2 : 0.8 : 1.0). Similarly, 
for 2 times elongations, the six asymmetric renditions include shapes for the a-axis (2.0 : 0.8 : 1.2 
and 2.0 : 1.2 : 0.8), b-axis (0.8 : 2.0 : 1.2 and 1.2 : 2.0 : 0.8), and c-axis (0.8 : 1.2 : 2.0 and 
1.2 : 0.8 : 2.0). 

Figures [4| [sj and |6] show the a-, b-, and c-axis elongations/reductions, respectively. Each figure 
shows X = 1.2 — 3.0 elongations {top spectra), the fiducial asymmetric bricks {middle spectra), and 
X = 0.8 — 0.1 reductions {bottom spectra). The fiducial asymmetric bricks {middle) are analogous 
to the fiducial cube in Figs. [2] and |3] {middle). For each elongation/reduction of a specific crys- 
tallographic axis, solid lines and dashed lines represent the two different asymmetric renditions 
(0.8 : 1.2 and 1.2 : 0.8, pertaining to the two unaltered axes). Certain spectral features exhibit 
distinguishable differences between the two asymmetric renditions and we denote these features as 
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being 'sensitive to the effects of asymmetry'. Other spectral features are indistinguishable for the 
two asymmetric renditions and are deemed 'insensitive to the effects of asymmetry'. Comparing 
and contrasting the asymmetric brick elongations/reductions in Figs. |4| [5| and [g] with the sym- 
metric cube elongations/reductions in Figs. [2] and |3j shows: 1) the trends in feature shape, peak 
location and relative strength that exist for the elongations/reductions of the cube also hold for 
elongations/reductions of both asymmetric renditions of the brick; and 2) that certain, but not all, 
spectral features of a given elongation/reduction of a specific axis exhibit differences in their feature 
shape, peak location, and relative strength between the two asymmetric renditions (0.8 : 1.2 versus 
1.2 : 0.8 for the axes that are not elongated/reduced). 

For the a-axis elongations of the asymmetric brick (Fig. [4| top), the 10.6 //m feature in the 
symmetric case (of the cube. Fig. [2]) 'splits' (bifurcates) into a 10.4 /im feature for the longer c-axis 
asymmetric rendition (0.8b : 1.2c) and a 10.6 ^m feature for longer b-axis asymmetric rendition 
(1.2b : 0.8c); the 11 fim feature is insensitive to asymmetry, i.e., the dashed and solid lines plot over 
one another. Furthermore, with a-axis elongations, the 11 fim feature shifts to longer wavelengths 
and the 10.4 or 10.6 nm features only slightly shift to shorter wavelengths and slightly decrease in 
relative strength, mirroring the cube-elongations (Fig. [2| top). 

For the b-axis elongations of the asymmetric brick (Fig. [5| top), the 9.7 fim feature in the 
symmetric case splits into 9.6 fim for longer a-axis asymmetric rendition and 9.8 ^m for the longer 
c-axis asymmetric rendition; the 11 ^m feature is insensitive to asymmetry. 

For the c-axis elongations of the asymmetric brick (Fig. [6| top), the features shortward of 
11 f^m show more complicated behavior than the a- or b-axis elongations and are divergent for the 
two asymmetric renditions (0.8a : 1.2b and 1.2a : 0.8b): for 'extreme' elongations, longer a-axis 
renditions (1.2a : 0.8b : Xc, dotted lines) develop 10.5 //m peaks whereas longer b-axis renditions 
(0.8a : 1.2b : Xc, solid lines) maintain a trough. For c-axis elongations, the 11 ^m feature is sensitive 
to asymmetry. Despite the variation in sensitivity to asymmetric renditions, all of the trends in 
spectral feature peak location shifts are equivalent to the symmetric case and the distinguishing 



features associated with elongations/reductions are present (§ 4.1.1, last paragraph). 



The spectral characteristics of the elongations with respect to the two asymmetric renditions in 
Figs. |4| [sj and |6] can be characterized by their primary axis elongation/reduction spectral feature 
characteristics coupled with spectral characteristics related to elongations of the second longest 
(1.2) axis. As already mentioned, for the X > 2.6 c-axis elongations in Fig. |6| the a- (dashed) 
and b-axis (solid) asymmetric renditions (1.2a : 0.8b and 0.8a : 1.2b) exhibit significantly different 
absorption characteristics between the 10 and 11 /xm features. The longer a-axis rendition of the 
c-elongated asymmetric brick has a distinct feature at 10.5 /im and a trough between the 10.5 /im 
feature and the 11 /um feature, which is similar to the a-axis cube-elongation's distinct 10.5 iim 
peak. The longer b-axis rendition of the c-elongated asymmetric brick has a broad feature centered 
near 10.1 fim and a shelf centered near 10.7 /xm that connects with the 11 iim feature, which is 
similar to the b-axis moderate cube-elongation's shelf feature near 10.8 /xm. 
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Fig. 4. — The 8 ~ 13 /^m (left panel) and 13 - 40 iim {right panel) absorption efficiencies normalized 
to the 11 fiui peak for the a-axis elongations {top curves) and reductions {bottom curves) for 
tri-axially asymmetric bricks compared to 'fiducial' bricks with axial ratios 1.0 : 1.2 : 0.8 and 
1.0 : 0.8 : 1.2 {middle curves). The dashed curves have the b-axis as the second longest (1.2) axis, 
and the solid curves have the c-axis as the second longest axis. The elongations are 1.6, 2.0, 2.6, 
and 3.0 times the 1.0 side of the 'fiducial' bricks, and the reductions are 0.6, 0.4, and 0.1 times the 
1.0 side of the 'fiducial' bricks. The 13 -40 fj,m Qabs values are multiplied by a factor of 1.5 for 
visual clarity. 
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Fig. 5. — The 8 ~ 13 (left panel) and 13 - 40 //m (right panel) absorption efficiencies normalized 
to tlie 11 /xm peak for the b-axis elongations (top curves) and reductions (bottom curves) for 
tri-axially asymmetric bricks compared to 'fiducial' bricks with axial ratios 0.8 : 1.0 : 1.2 and 
1.2 : 1.0 : 0.8 (middle curves). The dashed curves have the a-axis as the second longest (1.2) axis, 
and the solid curves have the c-axis as the second longest (1.2) axis. The elongations are 1.6, 2.0, 
2.6, and 3.0 times the 1.0 side of the 'fiducial' bricks, and the reductions are 0.6, 0.4, and 0.1 times 
the 1.0 side of the 'fiducial' bricks. The 13 -40 fj,m Qabs values are multiplied by a factor of 1.5 for 
visual clarity. 
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Fig. 6. — The 8-13 (left panel) and 13 - 40 //m (right panel) absorption efficiencies normalized 
to the 11 fj,m peak for the c-axis elongations (top curves) and reductions (bottom curves) for 
tri-axiahy asymmetric bricks compared to 'fiducial' bricks with axial ratios 0.8 : 1.2 : 1.0 and 
1.2 : 0.8 : 1.0 (middle curves). The dashed curves have the a-axis as the second longest (1.2) axis, 
and the solid curves have the b-axis as the second longest (1.2) axis. The elongations are 1.6, 2.0, 
2.6, and 3.0 times the 1.0 side of the 'fiducial' bricks, and the reductions are 0.6, 0.4, and 0.1 times 
the 1.0 side of the 'fiducial' bricks. The 13 -40 fj,m Qabs values are multiplied by a factor of 1.5 for 
visual clarity. 
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Feature sensitivity or insensitivity to asymmetry depends upon the specific axis that is being 
elongated/reduced. The variabihty in sensitivity can be understood through comparing symmetric 
elongations to the asymmetric elongations, as is demonstrated with the following two examples. 
The 11 fim feature is sensitive to asymmetry for a- and b-axis elongations, but not for c-axis 
elongation. Comparison to the elongations in the symmetric case indicates that the c-axis does 
not play a dominant role in determining the wavelength position of the 11 ^m feature, since in the 
symmetric case c-axis elongation only shifts the 11 /xm feature to shorter wavelengths by 0.1 nm 
over all elongations. For a- and b-axis elongations, however, the 11 fiin feature shifts by w 0.4 ^um 
in the symmetric case, and hence whether the a- or b-axis is longer in the asymmetric case plays a 
strong role in determining where the 11 fim feature peak is located. As with the 11 fim feature, the 
16 fim feature is sensitive to asymmetry for a- and b-axis elongation (Figs. |4] and [sj solid vs. dashed 
lines), but not for c-axis elongation. The symmetric elongations show that the relative strength of 
the 16 fim feature only significantly changes for c-axis elongations. Hence, the relative strength of 
the 16 /xm feature for a- and b-axis elongations depends on whether the c-axis is the second longest 
axis. The difference between this example and the previous one is that here the c-axis is important 
for the relative strength of the 16 fim feature, while for the previous example the c-axis is relatively 
unimportant in determining the position of the 11 fxm feature. A general rule is if a certain spectral 
resonance is sensitive to only a-, b-, or c-axis elongation/reduction and not the other two, then that 
feature will be sensitive to asymmetry for renditions that include it as a second longest axis. If the 
spectral resonance is sensitive to multiple axial elongations/reductions, then the feature will most 
likely be sensitive to asymmetry in all axial elongation/reduction cases, but to varying degree (e.g, 
the 10.5, 19, and 23.5 fiui features). 

Consider the fact that the spectral differences between the two asymmetric renditions for a 
specific crystallographic axis elongation/reduction do not systematically increase or decrease with 
increased elongation/reduction. For example, for a-axis elongations, the wavelengths of the split 
features do not change significantly with increased elongation for the (9.6, 9.8 ^m) and (10.7, 
10.5 /im) features for (b-,c-) second longest axis, respectively. The expectation is that if elonga- 
tion/reduction effects solely dominate the spectral features' shape, peak location, and size, that 
with increased elongation or reduction, the effects of asymmetry should lessen. This, however, is 
not the case, as the effects of the second longest axis are apparent in the most extreme elongations 
and reductions. The failure for the asymmetric effects to dampen out with increased elongations 
indicates that the elongations/reductions are the primary shape effects on the 8-40 //m spectral 
features and the effects of asymmetry must still be considered. Hence, the asymmetric triaxial 
shapes of forsterite crystals are important shape characteristics with respect to spectral feature's 
shape, peak position, and relative strength. DDA computations allow explorations of asymmetric 

. Since triaxial asymmetric 
shapes probably exist in lab samples and astronomical sources, the varying resonant features' sen- 
sitivities to the asymmetric renditions of the brick highlight the need to account for the asymmetric 
triaxial nature of forsterite in modeling and interpreting observational data. 



shapes or 'crystallographically anisotropic shape' (Takigawa et al.||2009 
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4-1.3. Prismatic and Complex Polyhedral Grain Shapes 

Grain shapes up to this point have been hmited to bricks (rectangular prisms), which have 
been employed to identify elongation/reduction and asymmetry of forsterite as important grain 
shape characteristics that affect the spectral features' shape, peak location, and relative strength. 
Brick shapes have each of their faces parallel to one of the crystallographic axes. Forsterite crystals 
in astrophysical environments, however, may not have their faces aligned in this manner. In a study 
of the crystalline silicate grains, including forsterite, found in Stardust tracks, the grain shapes are 
reported as euhedral, or rather, polyhedral shapes that have well-defined, sharp crystalline faces 
( Nakamura-Messenger et al.||2012 ). 

In this section, the DDA modeling exercises are extended to include non-brick polyhedral 
shapes including other prisms (triangular and hexagonal), pyramids (tetrahedra and rectangular 
pyramids or octahedra), and dipyramid^ (elongated triangular dipyramid, elongated hexagonal 
dipyramid, and dipyramidal brick). The same axial ratios are set across the grain shapes, which 
controls the effects of elongation/reduction and asymmetry. Hence, the observed spectral differences 
between sets of grain shapes with the same axial ratio are due to differences in the external shape. 

The 8-40 ;um spectral characteristics for sets of brick and non-brick shapes with similar 
crystallographic axial ratios are displayed in Fig. [7j For each set of shapes with similar axial 
ratios, all of the spectral features show remarkable similarities in feature shape, peak position, and 
relative strength, and each set is reflective of its elongation/reduction characteristics. Specifically, 
the central wavelengths, widths, and general asymmetries of each feature or feature complex are 
similar. Spectral differences, however, do exist, and they are primarily for the 10, 11, and 23 /xm 
features; from brick to non-brick the 11 /im feature shifts to longer wavelengths with the exception 
of the cube to hexagonal prism. Also, the shapes of the 23 ^m features vary. For the cube to 
octahedron and from the symmetric b-platelet to b-rectangular pyramid, the feature changes from 
trident to stair-stepped. For the c-axis elongated (c-long) brick to the c-axis triangular prism, the 
23 /im feature becomes stronger and peaks at a longer wavelength; note that the number of dipoles 
along the c-axis is greater for the triangular prism than for the elongated triangular prism because 
of the tips being incorporated into the c-axis elongation. Generally, the spectral differences between 
brick and non-brick shapes can be described as the non-brick crystalline resonances having weaker, 
rounder, and broader 10 and 11 ^va. features. These spectral differences suggest that external 
grain shapes with faces that are not solely parallel to the crystallographic axes are important 
shape characteristics. Moreover, these differences are typically greater between prismatic shapes 
and (di)pyramidal shapes suggesting that the presence of pyramidal structures is potentially an 
important shape characteristic. 

The effects of including pyramidal tips can be constrained through a comparison of the spec- 
tral features for prismatic shapes to their dipyramidal pairs of similar axial ratios, e.g., compar- 



^ Dipyramidal polyhedra are equivalent to bipyramidal polyhedra. 
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Fig. 7. — The 8-40 //m normalized to the 11 feature peak absorption efficiencies for rectangular 
prism (bricks) and non-brick shapes that have similar axial ratio. The axial ratios from top to 
bottom are approsimately a : b : c = (1 : 1 : 2), (1 : 1 : 1), (1 : 0.5 : 1), (2 : 1 : 1). Each set of curves 
contains a brick (solid curves) and a (di)pyramidal shape (dashed curves). All sets except the 
(1 : 0.5 : 1) also contain a non-brick prismatic shape. The 13 - 40 fim Qabs values are multiplied by 
a factor of 2 for visual clarity. The shape abbreviations are ETDP - Elongate Triangular Dipyramd; 
Tri. - Triangular; DiPyr - Dipyramidal; and Rect. - Rectangular; and Hex. - Hexagonal. 



ing triangular prisms to triangular dipyramids and to elongated triangular dipyramids (triangular 
dipyramids with trunks) with the same axial ratios. Such a comparison is shown in Fig. |8] for 
2 X elongated shapes for the a- (top), b- (middle), and c-axis (bottom) for the 8-40 spectral 
range. The prism - dipyramidal pairs presented in Fig [8] are the triangular prism - elongated trian- 
gular prism (ETDP), brick - dipyramidal brick (DPB), and hexagonal prism - elongated hexagonal 
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dipyramid (EHDPQ 

As with the brick to non-brick comparison, the 8-40 crystalhne resonance features for 
each prism - dipyramid pair are strikingly similar to one another with the characteristics of a-, 
b-, and c-elongation/reduction clearly identifiable. The positions and widths of the 10, 10.5, 16, 
and to a lesser degree, the 19 fiui features are the same, as expected from the elongation/reduction 
exercises on asymmetric bricks. The largest differences between the prism - dipyramid shapes are 
the peak location of the 11 feature, the width of the 11 and 23 fiui features, and the relative 
strengths of the 11, 19, and 23 fiui features. The dipyramidal shapes' spectral features generally 
have broader peaks and shoulders. Specifically, in all cases except the 2 x b-elongated brick - b- 
elongated dipyramidal brick (DPB) comparison, the 11 ^um feature is weaker in relative strength 
with the presence of pyramidal tips, and in all cases it is shifted to shorter wavelengths. The 
19 fj,m feature's behavior with the inclusion of pyramidal tips is dependent upon which of the three 
crystallographic axes is elongated. For example, the 19 //m feature systematically moves to shorter 
wavelengths and decreases in strength for the b-axis elongation case, while the 19 fim feature's 
secondary long-wavelength peak becomes a weaker shoulder for the c-axis elongation case. For a-, 
b-, and c-axis elongations, with the presence of pyramidal tips, the 23 /im feature always is at 
longer wavelengths for the pyramid compared to its dipyramidal pair. 

Comparing rectangular prisms to non-rectangular prisms and to dipyramids (Fig. [7| top, sec- 
ond, and bottom), the rectangular prisms (with faces parallel to crystallographic axes) have the 
23 //m feature terminating at the shortest wavelength. Inclusion of non-parallel faces extends 
the 23 ^m feature to longer wavelengths than for parallel faces, and to the longest wavelengths 
for dipyramids. While significant, this shape characteristic is not degenerate with the elonga- 
tion/reduction and asymmetry grain shape characteristics. Instead, it adds a layer of detail onto 
our understanding of what characteristics of grain shape determine the precise resonance feature 
shapes, peak positions, and relative strengths. The differences between the spectral features in 
Figs. [7| and [8] indicate that an important external shape characteristic is the presence of crystalline 
faces and edges that are not parallel to the crystallographic axes. 

4.2. Shape Classes: Columns, Platelets, and Equant 



The previous sections demonstrated that the effect of grain shape on the 8-40 ^m spectral 
features can be primarily described in terms of elongation/reduction with secondary spectral feature 
modifications due to asymmetry (non-elongated/non-reduced axes having unequal lengths) and 
non-parallel faces including tips on dipyramidal shapes. Specifically, the spectral features primarily 



comparison of a triangular dipyramid to elongated triangular prism [not shown] with the same axial ratio 
demonstrates that the two dipyramidal shapes are nearly spectrally indistinguishable from one another 
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are determined by the precise crystallographic axial ratio with the columnar, platelet, and equant 
grains' spectral shapes being clearly identifiable with non- linearly superimposed characteristics of 
the second longest axis and subtle, yet potentially significant, alterations to the features coming 



from non-parallel vs. parallel faces. In § 4.1.1 (last paragraph) we detailed a suite of diagnostic 
features for each of the a-, b-, and c-axis elongations. In both the symmetric and asymmetric 
cases, the behavior of the spectral features with elongation/reduction typically exhibited a distinct 
change for 'extreme' elongations {X > 2.6) and reductions {X < 0.2). Using these suites of 
diagnostic features, we define seven shape classes: a-, b-, and c-columns {X > 2.6 x); a-, b-, and 
c-platelets {X < 0.2 x); and equant (0.2 ^ X < 2.6). These shape classes serve as a basis set 
that can potentially be applied to astronomical observations to determine if the forsterite crystals 
are column-shaped ('whisker' or 'ribbon') or platelet-shaped. A practical application is using the 
shape classes to rule out 'extreme' elongations/reductions, i.e., reject 'whisker' or platelet shapes, 
and hence conclude equant shapes are present. 

The normalized absorptivities for the shape classes are shown in Fig[9j Each shape class is an 
average of the symmetric and asymmetric bricks within the defined elongation/reduction bounds. 
For each class, there are several 8-40 nm diagnostic feature shapes, positions, and relative strengths 
that make them distinguishable from one another. The a-columns have a 10.5 /^m feature that is 
equal to or stronger than the 10 fim feature, a birfurcated 23 //m feature, and an 11 fj,m feature 
that is at least twice as strong as all the other features. The b-columns have singular peaks for all 
the features, an increased absorptivity tail on the long wavelength side of the 11 fim feature, and a 
19 fim feature that is approximately twice as strong as the 23 /um feature. The b-columns are most 
easily identifiable by determining that the spectral features are not consistent with the other shape 
classes. The c-columns have a distinct 10-11 fxm feature complex that has a broad increased 
absorptivity shelf beginning near 10.0 ^m that is punctuated with distinct 10.5 and 11 //m peaks; 
a bifurcated 19 fj,m feature; and a singular 23 //m feature that is significantly enhanced in strength. 
The a-platelets have a broad 10.2 /xm feature coupled with a narrow 10.5 fim peak, and 19 and 
23 /im features that are nearly equal in strength and both enhanced in strength. The b-platelets 
have a broad 10.2 //m feature; lack a 10.5 /im peak; have a 19 //m feature with enhanced strength; 
and have a 23 /im feature with significantly enhanced strength. The c-platelets, of all the shape 
classes, have the largest separation between the 10 and 11 /xm peak positions with the 10 and 11 fim 
features at 9.5, 11.35 /im, respectively. The c-platelets have the 10 /im feature reduced in strength 
to the point of being a weak feature, and the 19 and 23 fim features are reduced and nearly equal 
in strength. 

The equant shape class in Fig. [9] is the average of all symmetric and asymmetric bricks with 
a-, b-, and c-axis elongations and reductions spanning 0.2 ^ X < 2.6 and with the non-altered 
axes having ratios of 0.8 : 1.2. Constituents of the equant shape class have crystallographic axes 
of nearly equal lengths so the features are sensitive to dominant effects of elongating/reducing the 
a-, b-, and c-axes and to the effects of the second longest axis. For example, the a-axis can be the 
second or the third longest axis: for c-elongated asymmetric bricks, the 10.5 iim feature is in the 



- 26 - 



I. 2 : 0.8 : 1.6 axial ratio brick but not in the 0.8 : 1.2 : 1.6 axial ratio brick, so the 10.5 fiui feature 
is present but is not dominant in the equant shape class. With the moderate elongations employed 
in the equant shape class, the 10 fim feature appears at 9.8 //m and the 11 //m feature appears at 

II. fj,m. The 19 //m feature is a singular peak. The 23 feature is a rounded version of the 
23 ^m trident feature, which is characteristic of the cube, as well as the other equant shapes of 
sphere and octahedron. The equant shape class spectral features in descending order of strength 
are: 11, 19, 23, 9.8, 16, 27, and 33 /zm. If spectral features are not attributable to any of the six 
extreme elongations and reductions (columns and platelets), then probably the equant shape class 
is relevant. Spectra of comets and protoplanetary disks are strikingly similar to the equant shape 
class (§ [5^ . 

The strength of the 10.5 um feature alone should not be used as diagnostic for the equant class. 
For equant shapes that are slightly elongated in the a-axis, either for the longest (Figs [2] and [3]) or 
second longest (Fig. |4]) axis, the 10.5 feature is enhanced. Even slight elongation {X = 1.2) of 
the a-axis produces a singular 10.5 fim peak. Alternately, an a-column identification can be made 
from a clear 10.5 ^m peak by examining the wavelength spacing between the 10 and 11 fiui peaks 
and the presence of a strong shoulder on the short wavelength end of bifurcated 23 ^m feature. The 
remainder of the diagnostic features of the equant shape class, however, remain clearly discernible. 
The equant shape class 23 fim feature is an excellent representation of the trident-shaped feature 
of the sphere, cube and octahedron triple-peak and a good generalization of the width and shape 
of the tetrahedron, albeit the tetrahedron is double humped. Generally, a-columns have a clearly 
bifurcated 23 fim feature while equants are not bifurcated. 



4.3. Effects of Grain Size on Resonant Features 



Grain size is an important factor when considering the spectral signature of forsterite grains 
in a natural setting. Within the Rayleigh domain, where the grains are small compared to the 
wavelength of light, the infrared features are strongly influenced by the resonances in the bulk 



refractive indices and are frequently referred to as 'surface modes' (Bohren & Huffman 1983: Min 



2009). As the wavelength of light approaches the grain size, or the 'size parameter' x = 2Tr x aeg/A 



approaches unity, there is a phase lag of the wave inside the grain because the wavelength of the 
radiation inside the grain is small; the wavelength in the grain is small because the refractive index 
at the strong spectral features can be very high. A feature becomes influenced by increasing grain 
size for increasing values of the product m x x = 27r x (m/A) (Min et al. 2005a; Min 2009), so 
shorter wavelength features (e.g., the 8-13 ^m features) are affected by increasing grain size at 
smaller grain sizes than the longer wavelength features (e.g., the 13 - 40 um features). A good 
example is the comparison between DDSCAT computations of ellipsoids and CDE, which is a time- 



independent solution applicable when the phase lag in the grain is negligible (Bohren & Huffman 



1983): Qabs agrees between DDSCAT and CDE for ellipsoidally shaped grains of sizes smaller than 
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fflcff — 0.5 fj,m, but for larger grains, the shapes of spectral features change with increasing grain 
size. Grain sizes in comets and laboratory samples typically span submicron to micron-sizes, so 
grain size effects are best explored with the DDA method. 

In this section, the effect of grain size on the 8-40 /im spectral features is investigated using 
DDSCAT for various grain shapes and sizes ranging from 0.1 - 3.0 /im in effective radius. In order 
to accomplish this goal, tetrahedra and elongated hexagonal dipyramids (EHDP) are computed 
with DDSCAT for ten grain sizes spanning effective radii Oefr = 0.1-1.0 with a linear step-size 



of 0.1 fim (Fig. 10). The absorption efficiencies are normalized to the 11 /im peak, the strongest 
peak for this range of grain sizes, to emphasize the change in the relative strengths of the features 
with changing grain size. The general trend in the 8-13 fj,m wavelength region is the peak positions 
of the spectral features shift to longer wavelengths as grain size increases, regardless of the external 
crystal shape. For increasing grain size there is also a decrease in the relative strength of the 10 /im 
feature, and in increase in relative strength for the 10.2 and 11.9 fim minor peaks. Understanding 
this behavior is essential to interpreting remote observations of comets since the 8-13 /im region of 
a comet's spectrum is the most commonly observed and utilized spectral region for diagnosing the 
silicate mineralogy ( Hanner et al.||1984 Wooden et al.||1999 ). In the 13 -!40 /im region, features do 



not exhibit significant shifts in peak wavelength position or feature shape, but do show a consistent 
increase in relative feature strength and broadening with increasing grain size. 



In Fig. 10, when the grain size increases from 0.1 to 1.0 /tm the relative strength doubles 
for all the 13 - 40 /xm features. As the grain size increases, the same trends in the 8-13 /xm 
wavelength region seen for the tetrahedral grains are evident for the EHDP grains. For both grain 
shapes, the locations of features at approximately 9.6 /im (shoulder for tetrahedra and peak for 
EHDP) and 11 /im shift to longer wavelengths. The amount of the shift of the 11 /im peak towards 
longer wavelengths is less, however, for EHDP (~ 0.1 /im) than for tetrahedra (~ 0.3 /im). The 
largest differences in feature shape with increasing grain size are that EHDP grains lose a short 
wavelength side shoulder on the 11 /im feature, and the largest change for the tetrahedra grains 
are the rounding out of the 10 iim peak and the 23 /im feature changing from a double-peak to a 
feature, which is better described as a short wavelength side shoulder with an increasingly strong 
long wavelength side peak. The difference in the amount of peak position shift between the two 
grain shapes indicates that grain size and grain shape effects are intertwined. The degree to which 
the peak positions shift with increasing grain size depends upon the shape characteristics of the 
grain. Nevertheless, the direction of the peak position shifts with increasing grain size remains the 
same regardless of specific grain shape. 

For forsterite, as grain size increases above ~0.4 /xm we see significant changes in the 8 - 
40 fim spectral shape due to the resonant features no longer being within the Rayleigh domain. 
Outside the Rayleigh domain, spectral features change their spectral feature shape, peak location, 
and relative strength depending on the wavelengths and the strength of the each of the resonances. 
Understanding these size dependency effects provides a diagnostic to determining the presence or 
absence of large grains in a sample of forsterite. 
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In the previous examples, the grain sizes were limited to a^s < 1.0 fj,m in size to clearly show 
the trends in feature location and strength with increasing grain size. Here, we extend the grain 
sizes to >1 fim to show the dramatic changes in feature shape that occur when well outside the 



Rayleigh domain. Figure 11 displays the grain shape of a forsterite dimension (FD) brick (a : b : c 
= 0.47 : 1 : 0.59, the axial ratio of forsterite's unit cell) for grain sizes 0.1 < a^s < 3.0 fim with 
a grain size resolution of Aa^s = 0.04 fim for 0.1 < a^s < 1.0 fim and Aa^s = 0.4 /im 
for 1.0 < Ocflf < 3.0 fim. The bottom spectra in Fig. [Tl] show all of the Qabs over-plotted 
and unnormalized to emphasize the changes to feature shape, peak location, and absolute feature 
strengths. The top spectra in Fig. 11 singles out Oefr = 0.1, 1.6, and 3.0 fim to more clearly elucidate 
large grain size spectral effects to the crystalline resonances. The un-normalized absorption spectra 
confirm the expectation that the shorter wavelength spectral features leave the Rayleigh domain 
at smaller grain sizes than the longer wavelength spectral features. Comparing the Oefr = 0.1, 1.6, 
and 3.0 nm isolated curves, shows that by Oefj = 1.6 /im, the 8-13 /xm features have significantly 
changed in their absorptive behavior, while the 13 - 40 fj,m features have primarily only increased in 
strength and remained similar in shape and wavelength location. At a grain size of Ceff = 3.0 /im, all 
but the 27 and 33 fim spectral features have left the Rayleigh domain, and the features have begun 
to significantly change in shape, peak location, and strength. In general, when a feature leaves the 
Rayleigh domain, the narrow crystalline resonances broaden significantly, become rounded, and 
decrease in peak absorption. Since the 11 /im feature diminishes in strength at smaller grain sizes 
than the 19 and 23 fim features, at intermediate grain sizes the strongest spectral feature changes 
from the 11 /xm peak to either the 19 or 23 ^m peak. Whether the 19 or 23 ^m peak is the 
strongest feature is dependent on the specific grain shape and grain size. In a power-law type grain 
size distribution, the smallest grains will dominate the surface area while the larger grains will 
dominate the mass. Despite the predominance by number of small grains, the presence of grains 
larger than ~1 /xm can be discerned in a grain size distribution. In astrophysical objects containing 
forsterite, e.g., the comae of comets, whether the 11 or 19/23 ^m feature is the strongest feature 
in the crystalline spectral component can serve as an observational spectral diagnostic to indicate 
whether or not a significant number of larger (>1 /xm) grains are present. 



4.4. Ensembles of Grain Shapes or Sizes 

When compared to the absorption characteristics of laboratory forsterite samples, our small, 
Ogfj = 0.1 fim grains have much narrower spectral features. The discrepancy is due to multiple 
factors, some of which are inherent differences between physical samples and our discrete models of 
defined shape and size target grains, and others of which are inherent to the nature of laboratory 
samples. In contrast to the single-size and single-shape DDA results we have presented thus far, 
any physical sample of forsterite contains both a distribution of grain shapes and a distribution 
of grain sizes. Additionally, laboratory prepared samples have grain properties that are related to 
the sample preparation techniques and the sample measurement techniques. Laboratory samples 
are ground down by either hand or ball-grinding techniques, which produce characteristic grain 
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sizes and shapes, and induce crystalline lattice distortions, all of which have been demonstrated to 
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Imai 


Imai et al. 



2009 ) that are not necessarily representative of astrophysical samples of forsterite. The difficulty in 
comparing model and laboratory results is further compounded by the measurement techniques. In 
order to measure the mass absorption coefficients for the ground samples, the laboratory samples 
either have to dispersed in an aerosol ( Tamanai et al.|[2006a[ 2009) or embedded in a 'transparent' 



medium of Cesium Iodide (Csl), Potassium Bromide (KBr), or polyethylene (PE) (Koike et al. 2003 



2006, 2010). The aerosol technique limits the grain sizes (oefr < 1.0 /um) as they pass through an 
impactor that attempts to separate agglomerated grains. Also, in the aerosol measurement chamber, 
the grains likely are electrostatically clumped. Both of these factors have been shown to affect 



the shape, peak position, and relative strength of spectral features (Tamanai et al. 2006a, 2009). 
Further, the embedding a sample in a medium for measurement has been shown to systematically 
shift the resonant features to longer wavelengths but where the amount of the wavelength shift 
varies depending on the feature ( Tamanai et al.||2006a ). 



Spectra of laboratory samples have broader and smoother spectral features than any single 
DDSCAT target grain of Oefr = 0.1. The spectral feature differences between laboratory samples 
and DDA modeled grains significantly lessens when we consider either an average of many brick 



shapes with different axial ratios (or shape ensemble. Fig. 12) or a single shape averaged and 
weighted over a grain size distribution (GSD) (or size ensemble. Fig. 13). The laboratory data 
that is lofted in air is the Alfa Irregular (Alfa Irr.) from Tamanai et al. (2006a), and the grain 
shapes in the sample are described as 'irregular'. The DDA modeled grains are an ensemble of 
bricks composed of all aes = 0.1 /xm grains, and are averaged over elongations/reductions without 
preference to any crystallographic axis. The two ensembles considered are: 1) an average of all 
the symmetric and asymmetric bricks ('all bricks ensemble') presented in § [4.1.1 and 4.1.1 and 
2) an average of all symmetric and asymmetric bricks where 0.4 < X < 2.0 ('shape restricted 
ensemble'). The 'shape restricted ensemble' is a brick-selected representation of the equant shape 
class described in § |4.2[ The two ensembles are compared to a single a-axis reduced asymmetric 
brick (0.4 : 1.2 : 0.8) for reference. Both the 'all bricks ensemble' and the 'shape restricted ensemble' 
have broader spectral features than the single brick shape. This demonstrates that when a variety 
of grain shapes are considered, the spectral features are broader and more comparable to physical 
samples of forsterite. 

The width of the features for the shape ensembles are still not as wide as the laboratory 
samples, but this exercise only considers grain shapes and restricts grain size to aeS = 0.1 /im, 



which is a much smaller grain size than the average grains in the Tamanai et al. (2006a) sample. 
When an ensemble of grain sizes or grain size distribution (GSD) is considered, the spectral features 
significantly broaden. Fig. 13 displays the Tamanai et al. ( 2006a[ ) Alfa Irregular laboratory sample 
data compared to a Manner GSD (Banner et al. 1984) with minimum grain size, ao 
peak grain size, a. 



0.4 lira; maximum grain size, a„ 



0.1 /im; 

1.0 //m; and power- law slope, N = 
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3.0. All of the absorption efficiencies for the GSD are DDA computed absorption efficiencies for 
tetrahedrally-shaped grains. These GSD parameters were chosen because they provide a realistic 



comparison to the Tamanai et al. (2006a) laboratory sample, which has unquantified grain size 
characteristics other than a maximum size limited by the aerosol injector at a^g = 1.0 /um. For 
comparison, the 'shape restricted ensemble' is included in Fig. [Tsjto exhibit the differences between 
the shape and size ensemble approaches. With the GSD approach, the spectral features are broader 
than for the shape ensemble. In general, the GSD spectral features are more comparable to the 
laboratory sample. The GSD approach also is able to better match the relative strengths of the 
laboratory spectral features than the shape ensemble approach for the single grain size of 0.1 /im, 
as is expected from the results of our study on grain size effects § |4.3[ 

Neither the shape ensemble nor the GSD approach replicate the precise features of the lab- 
oratory features, but both approaches demonstrate that including the realistic characteristics of 
a variety of grain shapes and grain sizes is able to broaden the features and make the relative 
strengths of the features more comparable to what is observed in laboratory experiments. It may 
be possible to provide high quality DDA-based model fits to laboratory data, but that is outside the 
goals of this investigation. Further, the approaches presented here only account for grain shape and 
grain size and do not attempt to account for the other previously mentioned grain characteristics 
including grain agglomeration in aerosol samples that also have an effect on the spectral features. 
Such an extended study probably would make the DDA-based models even more comparable to 
the laboratory data. 



5. DISCUSSION 

5.1. Discrete Grains: Sensitivity to Grain Shapes with DDSCAT 

We summarize each of the spectral resonances sensitivity to grain shape characteristics (elon- 
gation/reduction of the longest crystallographic axis, asymmetry, and complex faces) in Table [l] 
and note that some features are more sensitive to certain shape characteristics than other features. 
Certain spectral features exhibit varying degrees of overall sensitivity to grain shape such that some 
features are generally highly sensitive to all grain shape effects (10, 11, 19, 23, and 27 iJ,m features) 
or moderately sensitive (10.5, 16, 33.5 fj,m features). No single feature is weakly-to-not sensitive to 
all investigated grain shape characteristics, but certain features are selectively sensitive to certain 
grain shape characteristics. 

The feature sensitivity designations are strongly, moderately, and weakly-to-not sensitive to 
a particular shape characteristic, where the boundary between moderately sensitive and strongly 
sensitive is defined by whether or not the shape characteristic changes the spectral feature's shape 
as well as peak wavelength location and strength. Table [T] indicates that the 23 fj,m feature is the 
only feature that is sensitive to all of the grain shape characteristics, accentuating this feature's 
importance in analysis of thermal IR spectra of astrophysical objects containing forsterite. The 
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ll fim feature, on the other hand, only shows strong sensitivity to the presence of (di)pyramidal 
tips. The peak location of the 11 fj,m feature, however, varies significantly with all of the shape 
characteristics, but it is only with the addition of dipyramidal tips that the 11 fj,m feature exhibits 
a significant change in feature shape away from a well-defined singular peak. 

In addition to cataloging the general sensitivities of the 8-40 fim spectral features, we select 
a sample of the DDA grain shapes that represent the grain shape characteristics detailed in § [4] 



(shown in Fig. 14). For this sample, we identify the wavelength location of the peaks and shoulders 
of the spectral features (summarized in Table [2]) . The selected shapes are three columnar elon- 
gated triangular dipyramids (ETDPs) that are elongated in the a-, b-, and c-axes; three symmetric 
platelets that are reduced in the a-, b-, and c-axes; three equant asymmetric bricks [1.6 : 0.8 : 1.2], 
[0.8 : 1.6 : 1.2], and [1.2 : 0.8 : 1.6]; and four approximately 1:1:1 equant shapes: a tetra- 
hedron, octahedron, cube, and sphere. Such a selection includes elongated/reduced, asymmetric, 
and dipyramidal shapes, and a comparison to Fig. [9] demonstrates that each shape is well repre- 



sented by its shape class. Table [2] provides the quantitative basis and Fig. 14 provide a graphical 



representation to evaluate which spectral features are most sensitive to the effects of grain shape. 

Wavelength identification of the peaks and shoulders using an empirical algorithm (see Ap- 
pendix) gives a numerical description of the features shape through the location of the peaks and 
shoulders and quantitatively highlights the degree of variation in feature location caused by grain 
shape. The variation in the presence/absence of shoulders and minor peaks significantly varies over 
the presented grain shapes, highlighting the importance of grain shape in understanding the pre- 
cise spectral structure of observations. Without a proper treatment of the shape characteristics of 
forsterite in observations of comet comae and other astrophysical objects with crystalline forsterite 
spectral features, there is potentially a large problem in determining the precise mineralogy. Since 
forsterite is the strongest crystalline silicate feature in comet comae, inaccurate modeling of the 
forsterite spectral features potentially may lead to the necessity of additional minerals to account 
for the spectral differences between the forsterite model and the data. 



5.1.1. Summarizing Resonant Feature Sensitivity to Grain Shape 



Table [2] provides the quantitative basis and Figs. 14 provide a graphical representation to 



evaluate which spectral features are most sensitive to the effects of grain shape. This section 
analyzes the sensitivity of spectral features to grain shape for grains of a^s = 0.1 /um, with the results 
grouped according to the feature wavelengths. As a compliment to this section, refer to Table [T] 
that summarizes the sensitivity of each spectral feature to the effects of the shape characteristics 
of elongation, reduction, asymmetry, and the incorporation of pyramidal tips. 

The 10 - 12 jjLm complex has three strong features (10, 10.5, and 11 /im) and two weak features 
(10.2 ;um and 11.9 fim). The strong features exhibit a high degree of variation over the different 
grain shapes. For shapes where the c-axis is a strong component (c-axis columns and a-/b-axis 
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reduced platelets), the 10 /xm feature becomes blended with the minor 10.2 /xm peak and appears 

as a well-defined peak at ~ 10.2 /im. For grain shapes that have a triangular component (e.g., 
tctrahcdra, triangular dipyramid, and elongated triangular dipyramid), the 10 ;uni feature appears 
more as a broad shelf. This is a ubiquitous effect over all shapes with a 'triangular' shape, or rather, 
the shapes that have many acute angles between their edges and faces. 

The presence of the 10.5 //m feature is dependent on predominantly the a-axis, but also has 
a weaker manifestation due to the c-axis. When the a-axis is the longest or second longest axis, a 
10.5 /xm feature will be present. However, when the a-axis is minimized, the 10.5 /xm feature may 
not be present at all. With the a-axis ETDP, the presence of dipyramidal tips affects the shape 
of the 10.5 ^m feature. In this case, the broadening of the 11 /um feature due to pyramidal tips is 
extreme, and the 11 /xm feature connects with the 10.5 /xm feature turning the 10.5 /txm feature in 
to a shoulder. 

The 11 /im feature is a singular peak for nearly all of the grain shapes explored. The exception 
is the 'triangular' shapes, which have much broader 11 /xm features that are asymmetric primarily 
due to shoulders on the short wavelength side of the peak. The peak position of the 11 /xm feature, 
however, is very sensitive to grain elongation/reduction of a crystallographic axis. Over the shapes 
presented in this study, the 11 /xm feature peaks anywhere between 10.9 and 11.35 /tm. The 
elongation of the a- and b-axes is responsible for shifting the 11 /tm feature to longer wavelengths, 
while elongation of the c-axis is responsible for moving the 11 /tm feature to shorter wavelengths. 

The 16 jim peak is similar in shape across all of our grain shapes, but exhibits variation in 
peak position and strength due to c-axis contributions. 

The 19 iim feature exhibits significant variation in feature shape, peak positions, and strength, 
and is strongly sensitive to the elongation/reduction and asymmetry shape characteristics. The 
strong sensitivity to these two shape characteristics implies that the exact shape, peak locations, and 
strength are determined by the crystallographic axis ratio. Regardless of the exact crystallographic 
axis ratio, for columns and platelets, the 19 /tm feature maintains the general characteristics of the 
associated shape class. 

The 23 jim feature is the most sensitive to grain shape. It exhibits significant variation in 
feature shape, peak positions, and strength, and is strongly sensitive to the elongation/reduction, 
asymmetry, and non-crystallographically parallel faces shape characteristics. Over all shapes the 
23 /xm feature can appear as a single, double, or triple peak with complicated shoulder behavior 
depending on the grain shape characteristics. For highly symmetric 1:1:1 equant grains (i.e., 
sphere, cube, elongated hexagonal dipyramid), the 23 /tm feature appears at shorter wavelengths 
and tends toward appearing as a 'trident '-shaped triple peak. Breaking the high degree of symmetry 
for 1 : 1 : 1 equant shapes (e.g., tetrahedron), however, removes the 'trident '-shape, broadens, and 
shifts the 23 /xm feature to longer wavelengths. 

The weak 25 fim feature exhibits little variation over all the grain shapes. In cases when the 
c-axis is a strong component, or when the 23 /xm feature terminates at long wavelengths, the 25 /tm 
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feature will appear as a shoulder rather than a minor peak. 



The 27 fim feature exhibits only small variations in feature shape, and peak location for grain 
shapes that are not strongly influenced by the c-axis, i.e., c-columns or a-/b-platelets. In cases 
where the c-axis has a strong influence, the 27 ^m feature is strongly sensitive, and it appears as 
a distinctive 'shark tooth' shape that rapidly rises to the peak and slowly declines in absorptivity 
toward the 33.5 ^m feature. 

The 33.5 fj,m feature exhibits moderate sensitivity over the grain shape characteristics with 
the largest differences being in peak location and broadness of the peak. In Table [2j shapes with 
a notable influence from the a- and b-axes have a 33.5 fim feature that is indicated as having a 
fiat-topped two component peak. These 33.5 //m features also appear at longer wavelengths than 
features with strong c-axis influence. 

Out of all the 8-40 ^m forsterite spectral features, only the 16 and 33.5 /im features are 
moderately sensitive to all of the shape characteristics, implying that the character of forsterite 
grains' shape is a factor that cannot be overlooked when analyzing the other 8-40 fj,m spectral 
features. Interestingly, the 33.5 fim feature in laboratory data nearly ubiquitously exhibits two 



companion peaks or shoulders on the short and long wavelength sides of the main peak ( Tamanai 



et al.||2009t [Koike et al.||2010D . The observed 33.5 fim features for the comae of comets C/1995 01 



(Hale-Bopp) (Crovisier et al. 1997) at a heliocentric distance of 2.8 AU and 17P/Holmes during 
2007 outburst (Reach et al. 2010) also contain shoulders at ~32.8 and ~34.7 ^m. The presence 



of these additional shoulders on the 33.5 ^um features in physical samples and lack of them in our 
grain shape characteristic analysis suggests that these features are due to a grain characteristic 
other than those explored in this extensive investigation into forsterite shape characteristics. 



5.2. Model Applications 

5.2.1. Crystal Formation Environments 



Crystalline silicates in comets are the product of either gas phase condensation or annealing 
of amorphous grains. Condensation or annealing could have occurred in the inner 1 AU of the 



disk midplane where temperatures were > 1000 K (Harker et al. 2002 Wooden et al. 2005 2007 



Hanner &: Zolensky||2010 ). During the early high mass accretion phase (> 10 ^M0/yr) of the disk, 
the region hot enough to condense or anneal crystalline silicates could have extended out to 3 - 4 
AU dChick fc Cassen||1997t |Boss||1998[ [Beh et al.||2000D . There is still debate as to which of the two 
crystal forming processes was primarily responsible for the forsterite that is observed in the comae 
of comets. The grain shape of the crystals, however, has the potential to be a discerning factor in 
determining if the forsterite formed via condensation or annealing. 



Laboratory experiments that evaluate forsterite grain shape as the result of condensation and 
partial evaporation indicate that crystal shape and environmental conditions (temperature and 
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super-saturation) where the crystals form are causahy hnked (Tsuchiyama 1998; Bradley et al 



1983 Kobatake et al. 2008 Takigawa et al. 2009). One such laboratory condensation experiment 



is the Kobatake et al. (2008) experiment, where the investigators find that crystals that form from 



a rapidly cooled highly supersaturated silicate vapor on to metal substrates are characterized by 
'bulky', 'platy', 'columnar /needle', and droplet shapes for values of temperature and supersatu- 
ration, T and cr, of 1000-1450°C and < 97, 700-1000°C and 97-161, 580-820°C and 131-230, 



and <500°C and > 230, respectively (Kobatake et al. 2008). Their experimental columnar /needle 
shapes, which form by vapor-liquid-solid growth process at lower temperatures (<820°C), are ex- 
tended stacks of plates, where the extension is not correlated with an axial direction. Colum- 
nar/needles may be extended in the c-axis or a-axis direction, and can change directions, switching 
or combining a- and c-axis extensions, and appear off-kilter or a bit askew in shape. Other labo- 
ratory experiments find similar shape results, but there is not an exact agreement between crystal 
shape and environmental temperature or other conditions ( Tsuchiyama|1998 Takigawa et al.|2009[ ). 
Note that a preliminary report on condensation on to a forsterite single crystal substrate reveals 



rough spheroidal shapes (Tsuchiyama 2012), so the substrate composition and growth rate may 
affect crystal shape. 

While there is currently disagreement within the laboratory community between the forsterite 
crystal shapes produced and the environmental conditions, as the laboratory experiments become 



better and begin to converge, the seven spectral shape classes (§ 4.2) potentially can be applied 



to remote sensing observational data to determine the general grain shape characteristics present. 
Then, using the laboratory experiment results we can connect those shape characteristics to crystal 
formation environmental conditions. Such a link between spectral features, grain shape characteris- 
tics, and formation conditions is relevant to cometary comae, several of which have clear 8-40 ^m 
spectral features due to forsterite. Since cometary materials have remained unaltered since their 
formation and incorporation into comets, the determination of the crystal shape characteristics 
potentially provides a probe to the original crystal formation environmental conditions. 



As defined by Kobatake et al. (2008), the 'bulky', 'platy', and 'columnar/needle' shapes trans- 



late to the equant, platelet, and columnar shape classes (§ 4.2). Figure [9] shows the IR absorption 
efficiencies for the seven DDSCAT analog shapes that can be compared to the shapes of the con- 



densed grains in Kobatake et al. ( 2008 ) . Prom the descriptions and SEM images in Kobatake et al 



(2008), some of the bulky shapes are multi- faceted and others look like bricks with axial lengths 
similar to forsterite dimensions with a < b < c. We associate the bulky shapes, which are the 
highest temperature condensates (1270 - 1670 K), with the equant shape class. We associate the 
platy shapes, which are the second highest temperature (970 - 1270 K) condensates, with a-, b-, or 
c-platelets. We associate the a- and c-columnar shapes, which are coolest liquid-vapor-solid con- 
densates ( ), with a- and c-columns. B-columnar shapes do not appear to be one of the condensed 
shapes. A-platlets and a- and c-columns have the 10.5 ^m feature. A-platelets and c-columns 
have their 10 fim feature at 10.2 fim. A- and b-platelets and c-columns have a significantly en- 
hanced 23 //m feature (relative to other features). C-columns have their 11 /um feature shortward 
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of 11.0 /xm. Comparing the suite of diagnostic wavelength positions and relative strengths of the 
spectral features allows an assessment of whether equant, a- or c-columns, or platelets are prevalent 
in spectra of astrophysical sources. Identification of shape classes in astrophysical souces provides 
a potential means to probe the temperatures at which forsterite formed. 



5.2.2. Comets and Protoplanetary Disks 



Comet spectral energy distributions commonly show 9.8 and 11.2 ^um peaks without a strong 
and distinguishable 10.5 iim. peak ( Hanner &: Bradley||2004 > . In comet Hale-Bopp's SED observed 
near perihelion, the 10.5 /xm peak appears just as strong as the 11.2 /im peak but the crystalline 
contribution to the 10.5 /um peak is actually significantly weaker compared to the 11.2 /xm peak 
because the 10.5 /xm peak lies on the top of the amorphous silicate feature whereas the 11.2 fim 



peak lies on the trailing shoulder of the amorphous silicate feature (Wooden et al. 1999). Hence, 
the 10.5 fim peak in comets is not distinguishable and not strong, suggesting that the grains do 
not belong to the a-column, c-column, and a-platelet shape classes (§ 4.2 Fig. [9]). 



The cometary 11.2 /xm feature actually spans 11.05 - 11.2 ^m in high signal-to-noise spectra of 
comets Hale-Bopp ( [Crovisier et al.|1997[ [Wooden et al.|1999[ [Barker et al.|2002l|2004[ ), 17P/Holmes 



after its outburst (Reach et al. 2010), and 9P/Tempel 1 following Deep Impact (Lisse et al. 2006). 
The 11.05 - 11.2 /xm wavelength range for the peak excludes c-column grains as a component of 
Hale-Bopp's forsterite because c-column grains have their 11 /xm feature short-ward of 11.0 /im. 
These same three comets have their 10 /xm feature spanning 9.8 - 10.1 /xm. C-platelet grains have 
a very weak 10 /xm feature at approximately 9.6 /xm, so the comet grains likely do not belong to 
the c-platelet shape class. A significantly enhanced 23 /im feature relative to other features also 
is not observed for these comets so the grains likely are not c-columns or a-/b-reduced platelets. 
B-columns do not have a significantly enhanced 23 /xm feature and do have their 11 /im peaks near 
to 11.2 /xm, so b-columns are a likely spectral class in comets. The equant shape class produces 
the 11 and 23 /xm features at 11.0 and 23.5 /xm and the strongest feature is at 23.5 /xm. Specific 



axial ratios for equant shapes produce feature positions and relative strengths (Fig. 14) that are 
in good agreement with the spectra of these three comets: tetrahedra have their 11 /xm feature at 
11.05 /xm and a double-humped 23 /im feature with peaks at 22.85 and 23.40 /xm. Comparison of 
spectra of three bright comets with the seven shape classes suggests that the shapes of forsterite 
crystals in the comae of comets are equant or preferentially moderately elongated along the b-axi^ 
Exploration of the effects of crystal shape on spectral feature shape, peak positions, and relative 
strengths as well as exploration of the effects of grain size distributions sets the ground-work for 
fitting comet spectra with grain size distributions of equant or b-elongated forsterite grains and is 



within the scope of (Lindsay et al. 2013). 



^Forsterite dimensions is on the boundary between equant and b-columns. 
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A brief examination of external protoplanetary disk spectra show similar trends in spectral 
features as with comets, which historically has been shown through the classic comparison of comet 
Hale-Bopp with HD100546 ( [Malfait et ar||1998[ Fig. 4,). The 10.5 //m peak is not distinguishable, 
and the 23 /im peak is not overwhelmingly dominant (Olofsson et al. 2009; Oliveira et al. 2011). 
Likely, the grain shapes revealed by comet and disk observations are equant, and thus the precise 
spectral features are dependent on the crystallographic axis ratio. While this type of analysis cannot 
for certain determine the formation history of the crystalline grains, by the Kobatake et al. (2008) 
condensation experiment, the analysis suggests that the environmental conditions under which the 
crystals formed favored condensation at temperatures greater than 1000°C (1273 K). Condensation 
origin for forsterite crystals is strongly supported by laboratory examination of forsterite crystals in 
cometary IDPs ( Bradley et al.|[r983 Molster et al.|2003| >, which are described as having equant and 



tabular crystal habit (Bradley:2011, personal communication), and in Stardust samples (Brownlee 



et al.||2006| [Nakamura et al.||2008] [Tsuchiyama et al.|[2009l ). Hence, our shape classes can be 



used in an effort to ascertain the formation conditions of forsterite crystals that are observed in 
astronomical sources. 



6. SUMMARY AND CONCLUSIONS 



Through our DDA computations using the publicly available code DDSCAT v. 7.0 (Draine & 



Flatau 1994; Draine h Flatau 2008b), we identify three characteristics of grain shape that have 



distinct effects on the 8-40 //m spectral features of forsterite with respect to feature shape, peak 
position, and strength. In order of decreasing influence, the identified shape characteristics are: 
1) the extent to which a crystal is elongated/reduced along the longest/shortest crystallographic 



axis (§ 4.1.1); 2) the asymmetry of the crystal shape such that all three crystallographic axes are 



of different length, i.e. the crystallographic axial ratio (§ 4.1.2); and 3) the presence of crystalline 
faces that are not parallel to a crystallographic axis, e.g., (di)pyramids and non-rectangular prisms 



(§ 4.1.3). 



The first shape characteristic exercise, elongation/reduction, demonstrated that under 'ex- 
treme' elongation (> 2.6 x elongated along a single axis) and reduction (< 0.2 x reduced along 
a single crystallographic axis), the 8-40 ^m spectral features are distinct between the a-, b-, 
and c-axis elongations/reductions. The discerning spectral features are primarily the shape, po- 
sition, and strength of the 10, 10.5, 11, 19, 23, and 27 jim features. Using the distinctness of 
these spectral features between the elongation/reduction state, in § 4.2, we define seven spectral 
shape classes: 1-3) a-, b-, and c-axis columns; 4-6) a-, b-, and c-axis platelets; and 7) equant. 
Equant shapes are non-elongated/reduced to moderately elongated/reduced shapes that are nearly 
equal in all dimensions. Each of spectral shape classes has a suite of spectral features that can be 
used to potentially positively identify, or rule out, the elongation/reduction shape characterisitic 
of forsterite. 



The second shape characteristics exercise, asymmetry, demonstrates that the spectral char- 
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acteristics of the second longest crystallographic axis are important to the spectral signatures of 
forsterite. The spectral features for asymmetric bricks with three different axis lengths exhibit sig- 
nificant spectral differences between renditions of which crystallographic axis is the second longest 
axis. Regardless of how long or short the primary crystallographic axis is, spectral characteristics 
related to the second longest axis are non-linearly superimposed on the general spectral shape class 
template. The asymmetry effect particularly is apparent when the a-axis is the second longest axis. 
In this case, a 10.5 ^m feature that is characteristic of a-axis elongation is present. Such asymmetry 
effects can potentially lead to a misidentification of the spectral shape class, particularly for the 
equant shape class, which is the most sensitive to the effects of asymmetry due to the three axes 
being nearly equal in length. However, an examination of all of the 8-40 ^m spectral features can 
prevent a misidentification. 

The third shape characteristic, crystallographically non-parallel faces, demonstrates that the 
alignment of crystalline faces with respect to the crystallographic axis orientation can have sig- 
nificant spectral effects. This is accomplished by expanding the polyhedral shapes to include 
non-rectangular prisms and their (di)pyramidal pairs, which have the same axial ratios at the rect- 
angular prisms (bricks). The comparison between brick and non-brick shapes (Fig. [?]) confirm the 
expectation from the elongation/reduction and asymmetry results that shapes with similar axial 
ratios have similar spectral features. The agreement between spectral features is strong enough that 
regardless of the shape differences, the non-brick shape spectral features are clearly recognizable as 
belonging to their respective spectral shape classes. The main difference between the non-brick and 
brick shapes is the presence of crystal faces that are not parallel to the crystallographic axes imply- 
ing that the spectral differences are due to the presence of crystallographically non-parallel faces. 
The spectral effects of incorporating tips are subtle, and can be generally described as broadening 
the 10, 11, and 23 nm features. In summary, when the axial ratios are the same the differences 
between the spectral features are caused by the differences in grain shape. The effect due to the 
shape characteristic of crystallographic non-parallel faces also is strongest with shapes that have 
a 'triangular'-shape element, i.e., for tetrahedra, triangular dipyramids, and elongated triangular 
dipyramids. 

A quantitative comparison of the wavelength locations of spectral features for eleven DDSCAT 
polyhedral shapes and a sphere is given in Fig. 14 and in Table [2} Sensitivities of the spectral 



features near 10, 10.5, 11, 16, 19, 23.5, 25, 27, and 33.5 ^m to grain shape are summarized in Table[T] 



and discussed in § 5.1 When considering all shape characteristics, we recommend considering first 
the 10.5, 11 and 23.5 //m features, then the 10, 19, 25, and 27 /xm features, and lastly 16 and 
33.5 ^m features. 

A critical examination of the 8-40 /xm features with respect to grains size reveals that grain 
size alters the absorption resonances in distinguishable ways. In a grain size distribution of a specific 
grain shape, as a^s approaches 0.5 /um from below , the 8-40 /um spectral resonances begin to 
leave the Rayleigh domain. At aes > 1.6 fim, the changes in the shapes of the spectral features are 
no longer predictable based on the feature shapes of grains smaller than ~1.6 ^m. In general, as 
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grain size increases, the features diminish in strength and significantly change their shape starting 
at the 10-11 fj,m features and progressing to the longest wavelength features. 



The spectral shape classes defined in 



4.2| provide a practical means to analyze the shape 
Using the diagnostic spectral features 



characteristics of forsterite in comet comae observations, 
related to a-, b-, and c-axis columns/platelets, it is potentially possible to positively identify the 
shape characteristics of the forsterite in comet comae. By ruling out particular shape classes, the 
shape classes that possibly are present become apparent. In § |5.2.2[ we use the 10, 10.5, 11, and 
23 fim spectral features of those comet comae that have strong crystalline resonances to rule out 
the presence of platelet shapes and a- and c-axis columns. By the process of elimination, cometary 
forsterite likely is equant or b-columnar in shape. A more detailed analysis, however, is required 
to confirm this hypothesis. 



Based upon the results of condensation and partial evaporation experiments (Tsuchiyama 



1998t [Kobatake et"al]|2008t [Takigawa et"aI]|2009D , there is a causal link between the environmental 
conditions where forsterite forms and grain shape characteristics (§ 5.2.1). By comparison with 
condensation experiments by ( Kobatake et al.||2008 ), equant grains may be high temperature con- 
densates (1270 - 1670 K). Thus, knowledge of each of the grain shape characteristics can potentially 
provide insight into the formation environmental conditions of the Solar System's protoplanetery 
disk during the epochs of dust grain formation up to the time when these grains were incorporated 
into the small bodies of the Solar System. Parallel implications can also be drawn for external 
protoplanetary disks containing significant amounts of forsterite. 
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A. Algorithm for Identifying Peaks and Shoulders 

In order to consistently identify the peaks and shoulders in an absorption efficiency curve, 
we use an empirically derived algorithm that detects peaks and shoulders based on the curvature 
(second derivative) of the absorption efficiencies. The classifications of peak or shoulder are based 



upon the empirical Qabs curvature function, which is defined by Eqs. Al and 'K2 where 'sm' denotes 
the smooth function in the library of commands of the commonly used Interactive Data Language 
(IDL) (e.g., Albrecht et al.|[T997 ). If the strength of the curvature function is above a cut-off value 



for two or more consecutive wavelengths, then the related feature is identified by the algorithm 
and classified as a peak if there is a change in sign for the derivative of Qabs, and as a shoulder 
otherwise. The threshold value is chosen such that the shoulder of the 27 /im feature is included and 
less pronounced features are not identified. When identifying a peak, if a neighboring pixel Qabs 
value is within one percent of the maximum, the algorithm takes the average of the two wavelength 
positions as the peak location. 

There are cases when the curvature function has either a complex structure or is too narrow 
to return a proper identification of a peak or shoulder. A complex "M" shape in the curvature 
function occurs when a feature is broad. This is because the slope of the Qabs is rapidly changing 
as it approaches the peak or shoulder, slowly changing across the broad peak or shoulder, and 
then again rapidly changing on the far side of the peak or shoulder creating the m-structure in the 
curvature. For broad peaks, the algorithm selects the mid-wavelength position of the interval of the 
curvature function above the cut-off value for the feature in question. For the weaker shoulders, a 
linear fit is made between the endpoints of the range of the curvature function above the cut-off 
value for the feature in question. If the maximum of ratio of Qabs to this line-fit is greater than 2%, 
then the feature is identified as a weak shoulder. This empirical method produces identifications 
consistent with identifications by eye and yields a consistent set of wavelengths for the peaks and 
shoulders. 



Curvature Function = — q x irr-rriq} (Al) 



sm(sm(sm(sm({ 3), 20), 20), 20) 

Fig. [15] depicts how the peak identification algorithm is implemented for three grain shapes: 
sphere, tetrahedron, and dipyramidal brick. The DDSCAT Qabs values appear in black while the 
empirical curvature function is in red. In this figure, the curvature function has been normalized 
and plotted with a logarithmic ordinate to illustrate how the peaks and shoulders are identified. 
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The letters above the features indicate the peaks (P), broad peaks (B), shoulders (S), and broad 
shoulders (b). 
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Fig. 8. — The 8-40 absorption efficiencies for 2 x elongated prismatic shapes and their 
dipyramidal pair. The top panels are a-axis elongation, the middle panels are b-axis elongation, 
and the bottom panels are c-axis elongation. The prism - dipyramid paris are from top-to-bottom 
in each panel: emphblues triangular prism - elongated triangular dipyramid (ETDP); blue-greens 
brick - dipryamidal (DiPyr.) brick; and reds hexagonal prism - elongated hexagonal prism (EHDP). 
The 13 -40 //m Qabs values are multiplied by a factor of 2 for visual clarity. 
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Fig. 9. — The 8-40 normalized to the 11 //m peak absorption efficiencies for the seven spectral 
shape classes defined in text. The column shape classes are averages of all the symmetric and 
asymmetric bricks that are elongated by at least a factor of 2.6. The platelet shape classes are 
averages of all the symmetric and asymmetric bricks that are reduced by at least a factor of 0.2. 
The equant shape class is an average of all the symmetric and asymmetric bricks between elongation 
and reduction factors of 2.6 and 0.2. 
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Fig. 10. — The normalized to the 11 peak absorption efficiencies from 8-13 /im (left panels) 
and 13 - 40 nm (right panels) for DDA computed elongated hexagonal dipyramids (EHDPs; top 
panels) and tetrahedral (bottom panels) forsterite. The effective radius is varied from 0.1 /im to 
1.0 /xm over 10 grain sizes Unearly distributed with a 0.1 micron step-size. 
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Fig. 11. — The 8-40 //m absorption efficiencies for DDA computed forsterite dimension bricks 
(asymmetric bricks with a : b : c = 0.47 : 1 : 0.59) forsterite grains for grain sizes Oeff = 0.1 - 3.0 /im 
with a hnear step size of 0.04 //m for a^s = 0.1 - 1.0 /im and 0.4 ^m for = 1.0 - 3.0 //m (bottom 
spectra). The top spectra isolate the grain sizes, aes = 0.1, 1.6, and 3.0 //m to clearly demonstrate 
the effects of grain size as the spectral features leave the Rayleigh domain. 
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Fig. 12. — The normalized to the 11 fj,m peak absorption efficiencies for a single 0.4 : 1.2 : 0.8 
asymmetric brick (dashed green), the 'shape restricted ensemble' (solid red) of all symmetric and 
asymmetric bricks where 0.4 < X < 2.0, the 'all bricks ensemble' (solid blue) of all symmetric 



and asymmetric bricks, and the Tamanai et al. (2006a) aerosol technique measured Alfa Irregular 



laboratory sample (solid black). The widths of the features increases from the single shape to the 
'shape restricted ensemble' to the 'all bricks ensemble'. 



- 51 - 



O 

T3 
O 
_N 



- Tamanai et al., 2006a Alfa Irr. 

- 0.4 - 2.0 cubes and asym. bricks 
-Tetra GSD: ao=O.V(m, a =0.4^m 




9 10 11 12 
Wavelength (wm) 



13 15 20 25 30 
Wavelength (wm) 



35 



40 



Fig. 13. — The normalized to the 11 /.tm peak absorption efficiencies for a single 0.4 : 1.2 : 0.8 
asymmetric brick (dashed green), the 'shape restricted ensemble' (solid red) of all symmetric and 
asymmetric bricks where 0.4 < X < 2.0; a Hanner grain size distribution (GSD) average with 
parameters oq = 0.1 /xm, Op = 0.4 fim, Omax = 1-0 /^m, and power law slope, N = 3.0; and the 



Tamanai et al. (2006a) aerosol technique measured Alfa Irregular laboratory sample (solid black). 
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Table 1: Shape Characteristic Feature Sensitivities 



Characteristic 


10/im 


10.5/im 


11/im 


16/im 


19/im 


23/im 


27;um 


33.5/im 


a-elongation 


/ 


/ 


/ 




/ 


/ 


/ 


/ 


b-elongation 


/ 


X 


/ 




/ 


// 




/ 


c-elongation 


// 


// 


/ 


/ 


// 


// 


// 




a-reduction 


// 


// 


/ 


/ 


// 


// 




/ 



b-rcduction 


// 


// 


/ 


/ 


// 


// 


/ 


/ 


c-rcduction 


/ 


X 


/ 


/ 


// 


// 


/ 


/ 


elong./reduc. 


// 


// 


/ 


/ 


// 


// 


// 


/ 


Asymmetry 


/ 


// 


/ 


/ 


// 


// 


// 


/ 


Pyramidal Tips 


/ 


/ 


// 




/ 


// 


/ 


/ 



- Not-to-Weakly Sensitive: Little to no change in feature 

/ Moderately Sensitive: PeaJc location changes and minor strength changes,i-e., no shape change 
/ / Strongly Sensitive: Feature shape, location, and strength changes 
X: Feature not present 
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Fig. 14. — The normalized to the 11 /im peak absorption efficiencies for a grain shape characteris- 
tic representative sample of the DDSCAT computations presented in this paper. The grain shapes 
include three columnar shapes (ETDP^, ETDP^, and ETDPc; Elongated Triangular Dipyramids), 
three platelet shapes (Symmetric a-, b-, and c-Platelets) , three asymmetric equant shapes (asym- 
metric bricks), and three 1:1:1 equant shapes (tetrahedron, cube, and sphere). This set of grain 
shapes represents all seven shape classes and all three grain shape characteristics presented in this 
paper. A summary of the wavelength positions of the spectra's peaks and shoulders is in Table [2] 
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Table 2: Resonance Peak Locations 



Shape 


10/im 


10.5/um 


11/nm 


ii.g^m 




19;im 


23^m 


25/um 


27jUm 


33/im 


a-axis Elongated 

Triangular 

Dipyramid 


9.60 

9.80 


10.20 
10.50 


10.80 
11.25 


12.00 


15.60* 
16.00 


18.95 

19.90 


22.40 
23.90 


25.40 


27.45* 
27.80 

28.20* 


[33.60*, 
34.20+] 


b-axis Elongated 

Triangular 

Dipyramid 


9.60 

9.80 


10.20 


10.80* 
11.20 


11.85* 


15.60* 
16.00 


18.40^ 
19.35 


21.15 
23.35 
22.95 




27.40 


[33.55*, 
33.975] 


c-axis Elongated 

Triangular 

Dipyramid 




10.15 

10.50 


11.00 

11.50 


11.85 


16.40 


17.95* 
18.70 

19.55 


21.45 
23.35 


25.20 


26.60 

27.35* 


33.60 


Brick 

Symmetric 
a-platelet 




10.20 

10.50 


10.85 
11.35 




16.55 


17.80 
18.30 
19.675 
21.00 


21.80* 
23.40* 
23.85 


25.20 


26.625 

27.80 


34.25+ 


Brick 

Symmetric 
b-Platelet 




10.20t 

10.50 


11.35 


12.00 


16.55 


17.75 
18.35* 
19.00 
19.90 


22.20+ 
24.075 


25.60 
26.625 


27.80* 
28.225 


34.50 


Brick 

Symmetric 
c-Platelet 


9.40 


10.20 


11.35 


12.00 


15.55 


19.00 
19.60 

19.90* 


21.05+ 
22.175 
23.30 
24.20 


25.525 


28.225 


34.425 


Asymmetric 
Brick 

1.6 : 0.8 : 1.2 


9.80 


10.20 
10.50 


11.15 

11.50 


11.95 


16.05 


18.05 
18.75 

19.00* 


21.775 
22.55 
23.71 


25.00* 
25.325 


26.45 
27.45 

27.80* 


[33.60, 
34.00] 


Asymmetric 
Brick 

0.8 ; 1.6 : 1.2 


9.80 


10.20 


10.75 
11.05 

11.50 


11.85 
12.60* 


16.05 


18.10 
19.10 


21.15+ 
22.70 

.25.55 


25.15+ 


26.50* 
27.35 


[33.625, 

34.00] 


Asymmetric 
Brick 

1.2 : 0.8 : 1.6 


9.90 


10.20 
10.50 


11.025 

11.50 


11.90 


16.20 


18.05 
18.575 

19.30 


21.525 
22.80 
23.40 


25.00+ 
25.325 


26.55 
27.40 

27.80* 


[33.60, 

34.00] 


Tetrahedron 


9.70 
9.85* 


10.20 
10.50* 


10.90 
11.05 


11.95 
12.60* 


16.075 


18.75* 

19.00 

19.90 


21.80 

22.85 

23.40 


25.30 


26.55* 
27.475 


[33.60*, 
34.00] 


Cube 


9.80 


10.20 
10.65 


11.00 

11.50 


11.90 

12.60* 


16.00 


18.65 

19.05 

19.90 


21.475 
22.575 
23.40 


25.15+ 


26.45* 
27.375 


33.60 

34.00* 


Sphere 


9.80 


10.20 
10.70 


10.90 

11.50 


11.90 
12.10* 


16.00 


18.40 

19.00 
19.90 


21.40 

22.425 

23.35 


25.15^ 


26.45* 
27.35 


33.55 



Notation: major peak, peak, shoulder, * - minor shoulder, f - Broad PeaJc 
[ , ] - Flat-topped, two component peak. 
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Fig. 15. — The curvature function (red) compared to the DDSCAT Qabs values (black) for a 
tetrahedron (top), dipyramidal brick with forsterite dimension (middle), and sphere (bottom) grain 
shapes. The curvature function is defined as the square of the Qabs divided by the smooth of the 
Qabs multiphed by the negative of the second derivative of the Qabs- The dotted red line indicates 
the curvature strength cut-off value, such that features with curvature function value lower than 
the cut-off are not identified. The letters above the features indicate the peaks (P), broad peaks 
(B), and shoulders (S) 



